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Abstract 

In a spacetime with no global timelike Killing vector, we do not have a natural 
choice for the vacuum state of matter fields, leading to the ambiguity in defining the 
Feynman propagators. In this paper, choosing the vacuum state as the instantaneous 
ground state of the Hamiltonian at each moment, we develop a method for calculating 
wave functions associated with the vacuum and the corresponding in-in and in-out 
propagators. We apply this method to free scalar field theory in de Sitter space 
and obtain de Sitter invariant propagators in various coordinate patches. We show 
that the in-out propagator in the Poincare patch has a finite massless limit with the 
full de Sitter invariance preserved. We argue and numerically check that our in-out 
propagators coincide with those obtained by the path integral with the standard ie 
prescription. We also show that the in-out propagators satisfy Polyakov's composition 
principle. Several applications of our propagators are also discussed. 
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1. Introduction 



In a spacetime with no global timelike Killing vector, we do not have an established prescrip- 
tion to define the vacuum state of matter fields, along with which the propagator cannot be 
defined uniquely, de Sitter space is a typical example of such spacetime, and various vacua 
have been studied throughout the decades. Among them, the Euclidean vacuum (or the 
Bunch-Davies vacuum [I]) is de Sitter invariant and satisfies the Hadamard condition. This 
vacuum is regarded as the most natural vacuum in cosmology and is used to describe the 
physics in the inflationary era. Other de Sitter invariant vacua, called a- vacua [2j [3], have 
also been studied frequently (a = — oo corresponds to the Euclidean vacuum). However, for 
now, there is no conclusive reason to single out one vacuum. 

An interesting feature of de Sitter space is its thermodynamic property. As is pointed 
out in jl], a particle detector staying in de Sitter space behaves as if it is in a thermal 
bath with the temperature T = 1/2tt£, where I is the de Sitter radius. The calculation 
crucially depends on the choice of the vacuum at the initial time when the detector starts its 
operation. For example, the a- vacua do not predict the thermal behavior for a ^ — oo [5]. In 
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this sense, the choice of the vacuum is also important in understanding the thermodynamic 
behavior of curved spacetimes. 

In this paper, we define the vacuum of a free scalar field as the instantaneous ground state 
of the Hamiltonian at each moment. We develop a method to calculate transition amplitudes 
during finite time intervals for quantum mechanics with time-dependent Hamiltonian, and 
define the propagators as the limit of two-point functions when the initial and final times 
are sent to the past and future infinities^ In our method, wave functions associated with 
the vacuum are automatically determined with no need to consider asymptotic boundary 
conditions such as positive-energy conditions. 

We apply the method to construct various Feynman propagators in de Sitter space. We 
treat the principal (m > (d — l)/2) and the complementary series (m < (d — l)/2) at the 
same time, and show that the obtained in-in and in-out propagators always take de Sitter 
invariant forms. Furthermore, we show that our de Sitter invariant in-out propagator has 
a finite massless limit in the Poi.care patch! This is in contrast to the in-in propagators, 
for which the no-go theorem is known that there is no de Sitter invariant Fock vacuum for 
massless scalar fields [3]. 

We argue and numerically check that our in-out propagators coincide with the propaga- 
tors obtained by the path integral with the standard ie prescription. This result is consistent 
with the de Sitter invariance of the propagators since the corresponding path integral is per- 
formed over a patch which is invariant under the action of de Sitter group. Moreover, 
our in-out propagators are shown to satisfy the composition principle [9], which has been 
claimed by Polyakov recently as the property to be satisfied in order for the propagator to 
be interpreted as representing the sum over paths of a particle moving in a spacetime. 

This paper is organized as follows. In section [21 we develop a general framework to 
calculate the wave functions and propagators for quantum mechanics with time-dependent 
Hamiltonian. In section [3] we demonstrate how our prescription works in the simplest space- 
time, Minkowski space. Another well-studied example of asymptotically Minkowski space is 
investigated in Appendix [Bj We then analyze the de Sitter case in section |H After giving a 

1 There had been a study to define the vacuum as the instantaneous ground state, which is sometimes 
called the instantaneous Hamiltonian diagonalization method (see, e.g., [SI [71 [8] and references therein). Our 
framework has the same principle as that of the method in determining the vacuum, but has an advantage 
over the method in that it can determine an explicit form of various propagators for finite time intervals 
and can be applied for a wide class of nonstatic spacetimes. 

2 The in-out propagator in the global patch still diverges in the massless limit. 
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brief review on the geometry of de Sitter space in subsection I4.1[ we discuss the propagators 
in the Poincare patch in subsection 14.21 and those in the global patch in subsection 14.31 For 
the both patches, we first make a mode expansion of a scalar field and calculate the propa- 
gators for each mode. We then make the summation over modes to obtain the propagators 
in spacetime, both of the in-in and in-out types. The obtained propagators are found to be 
invariant under the full de Sitter group. In section [5] we argue and give an evidence based on 
numerical calculation that our in-out propagators coincide with those obtained by the path 
integral with the standard ie prescription. In section [6] we prove that our in-out propagators 
have the heat-kernel representation, which means that the propagators satisfy Polyakov's 
composition principle [10]. Section [7] is devoted to discussions and conclusion. We collect 
miscellaneous formulae in appendices. 

2. General framework 
2.1. Setup 

In this paper, we consider quantum theory of a free scalar field <j)(x) living in a d- dimensional 
curved spacetime with background metric g^J^ 

S[<Kx)\ = - \ J d d x yf^ (<T d^drf + m 2 2 ) . (2.1) 

We assume that the spacetime is globally hyperbolic and the foliation of spacetime (i.e., the 
set of timeslices) is already specified. We denote the temporal and spatial coordinates by t 
and x, respectively, and the spacetime coordinates by x = (a; M ) = (t,x) (/x = 0, 1, . . . , d— 1). 
We further assume that the metric has the following form: 

ds 2 = - N 2 (t) dt 2 + A 2 (t) hij(x) dx i dx j (i, j = 1, . . . , d - 1) . (2.2) 

The action is then written as 

5=1 J dt J d^xVhNA'*- 1 ^N~ 2 d t (Pd t ^ + A^ 2 <j)A d ^(P-m 2 (P 2 ^ , (2.3) 

where y/h = W det hij , and A^_i = (1/yh) di(\fhh 13 ' dj) is the Laplacian for the spatial 
metric ds 2 i _ l = hij(x) dx l dx J . We have neglected the surface term coming from integration 

3 The metric has the signature (—,+,...,+). 
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by parts. We introduce a complete system {Y n (x)} of orthonormal eigenfunctions of A^-i 
satisfying 

A d _i Y n (x) = - X n Y n (x) , J d d ~ l x y/h(xjY n {x) Y n ,{x) = 5 nn , , (2.4) 
and make the mode expansion of the scalar field as 

(j>{x) = <j)(t, x) = J2 Mt) Y n (x) . (2.5) 

n 

The action can then be written as the sum of the actions for mode functions {4> n (t)} '■ 

S = J2 [ <ttL n (<f> n {t),Mt),t) (2.6) 

with 

L n {(p n ,(j) n ,t) = — -0„ <j) n , (2.7) 



p (t) = N-\t)A d ~\t), u n (t)=N(t) ^\ n A- 2 {t)+m^. (2.8) 

This shows that the n th mode function 4> n (t) behaves as a quantum oscillator with time- 
dependent mass p(t) and frequency oo n (t), and is described by the following time-dependent 
Hamiltonian in the Schrodinger picture: 

H n , s (t) = H n {<p n ^ 7r n , s , t) = ^< g + P(t) if (t) (2.9) 

where the suffix s indicates that the operators are in the Schrodinger picture. Thus, the 
theory is reduced to quantum mechanics of a set of independent harmonic oscillators with 
time-dependent parameters: 

H.(t) = £ = Zy^<* + Pj ^ 1 <s] • (2-10) 

n n ^ 1 

The quantization is accomplished by setting the commutation relations 

[<fin,s, 7Tm,s] = i &n,m , [<Pn,s, 4>m,s] = = [lT n , s , 7T m)S ] . (2.11) 

In the following subsections, we develop a general theory to describe the time evolution of 
states for quantum mechanics with such a time-dependent Hamiltonian. 

We here make a comment on a subtlety existing in field redefinitions. By transforming 
the mode function <j> n {t) to Xn(t) = P~ 1/2 (t) <p n (t) = e a{t) (j) n {t) , one can make the coefficient 
of the kinetic term to unity: 

S[Xn(t)] = [ tf dt \ [xl(t) - Q 2 n (t) xl(t)} , (2.12) 
Ql(t) = ul{t) - (a{t)f - a(t) . (2.13) 
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However, it can often happen that fi^(i) takes a negative value for some region of m 2 even 
though the original co 2 (t) is strictly positive^ Although the physics should be the same 
for the two descriptions using <f) n (t) and Xn{t), the inverted harmonic potential for Xn(t) 
can easily cause a catastrophe when making an analysis based on an approximation such 



as the WKB approximation. In order to avoid this subtlety 
symmetry manifest), we will not make such transformations^ 



and also to keep the original 



2.2. Quantum mechanics with time-dependent Hamiltonian 

For brevity, in the following discussions, we omit the mode index n and denote the canonical 
variables {<p n ,s, ^n,s} in the Schrodinger picture simply by {q s ,p s }- Our Hamiltonian then 
takes the form 

H s (t) = H(q s ,p s ,t) = ^ f) P 2 s + ^^q 2 , (2.14) 
and the system is quantized by setting the commutation relation 

[q.,P.] = i. (2.15) 

As we will see below, in order to discuss the behavior of states near the temporal boundary 
in a well-defined manner, we need to let the frequency u(t) have a small negative imaginary 
part without breaking the symmetry existing in the background spacetime. This can be 
best realized by replacing m 2 in ( 12. ip by m 2 — ie (e > 0). Note that the Hamiltonian H s (t) 
is no longer hermitian when e ^ 0. 

We denote by T s the time at which quantization is carried out in the Schrodinger picture. 
The Hilbert state Ti = {\ip)} with the hermitian inner product (1^1,1^2) is then constructed 
on the timeslice at t = T s , and the dual space H* = is defined with respect to the 

hermitian inner product with the rule (ipi\ = (|^>i)) , i.e., (^lld^)) — (^i?^)- The time 
evolution of a state £ "H is governed by the Schrodinger equation 

d t \iP,t) = -iH s (t)\*l;,t) (2.16) 

4 A typical example is a scalar field in the Poincare patch of de Sitter space. One can easily see that 

f2^ (t) can be negative when the mass is small and <f> n (t) represents a mode of long wave length. 

5 The coefficient of the kinetic term can be set to unity also by making a transformation of the time 

coordinate. We will see that physical quantities do not change under the transformation (see the last 

paragraph of subsec. I2.4[) . 
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with the initial condition \ip,T s ) = The Schrodinger equation can be integrated to the 
form 

\1>,t) = U{t,T,)\iJ>), (2.17) 

where U(t,T s ) is the time-evolution operator expressed as the time-ordered exponential of 
H 8 (t), 

U(t,T s ) =Texp(-i J dfiZ,(f)) 

= lim (1 - iAt N H s (t N )) (1 - iA^-i# s (^-i)) • • • (1 - iAhHsih)) 

At k — >0 

(t = t N > tjv-i > ■ ■ ■ > h > t = T s ; At k = t k - t fc _i) . (2.18) 
The hermitian conjugate of \ip,t) is given by 

(,/>,t\ = (,l>\Ut{t,T a ). (2.19) 

In addition, we introduce a one-parameter family of states for a given state G %* as 



i^,t\ = ^\U-\t,T s ), (2.20) 

which satisfy 



d t (ij,t\ = +i(ij } t\H s (t) J &,T S \ = (V| . (2.21) 



Note that the paring of (^il and | ^2) does not change under the time evolution, 



<Vi,*|-0 2 ,t> = <-0il-02> , (2.22) 

although this is not the case for (ipi,t\ip2,t) when u(t) has a nonvanishing imaginary part 
because the time evolution operator is then not unitary, £7 _1 (t,T s ) ^ W(t,T s ) . 

The spectrum of the Hamiltonian H s (t) can be easily found by introducing, as usual, a 
pair of operators^] 



m ^,e^ q ^_l_ p ,. ( , 24) 



6 They are Schrodinger operators and the time dependence comes only through the parameters p(t) and 
w(f). 



We call a s (t) and a s (t) the annihilation and creation operators at time t. Note that a s (t) ^ 
a' s (t) when u(t) has a nonvanishing imaginary part. From the commutation relation (|2.15p . 
we have 

[a s {t),a s {t)} = 1. (2.25) 
The Hamiltonian ( I2.14p can then be rewritten as 

H s {t) = ^ [a s (t) a s (t) + a s (t) a s {t)} = u(t) [a s (t) a s (t) + i] . (2.26) 

We define the state \0 t ,t) as that which vanishes when acted on by a s (t): 

a 8 (t)\O t ,t) = 0. (2.27) 



Similarly, the state (0 t ,t| is defined a 



(0 u t\a s (t) = 0. (2.28) 



Then the left and right eigenstates of H s (t) are given by 

t (M[= 4 (Ml [a s (t)]\ \n,t) t = ^ [a s (t)TKt), (2.29) 

which satisfy 

t J^J\H s {t) = (n+-)w(i) t (M[, ^(t)|n,t) t = (n+ -) |n,t) t . (2.30) 

We call |0t, t) = |0, t) t the ground state (or the vacuum) at time t, since this is the minimum 
energy state when e = . 

It is important to note that the ground state at time t', \0 t i,t'), is generically different 
from the state \0 t ,t'); the latter is obtained as the time evolution of the ground state \0 t ,t) 
at time t, \0t,t') = U(t' ,t)\0t,t) (see Fig. [1]). Since the Hamiltonian is already specified at 
each time, the vacuum state is uniquely determined, and there is no freedom to introduce 
other vacuum states through Bogoliubov transformations without changing the timeslices 
(or Hamiltonian). 



2.3. Heisenberg picture 

Now we move from the Schrodinger picture to the Heisenberg picture. Given a Schrodinger 
operator O s (t) [possibly depending on t through the parameters involved when constructing 

7 Note that (0^| = (0 t ,T s | may differ from (0 t | = (0 t ,T s | since a s {T s ) ^ a^(T s ) when Imw(i) ^ . 
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Figure 1: Time evolution of states. The system is quantized in the Schrodinger 
picture on the timeslice at T s . \0 t ,t) is the state annihilated by the Schrodinger 
operator a s (t), while |0 f ) is the state annihilated by the Heisenberg operator a(t). 



the operator as in Eq. (12. 23ft ]. we define the corresponding Heisenberg operator as 

0(t) = U-%T a )0 9 (t)U(t t T 8 ), (2.31) 
which satisfies the Heisenberg equation, 

6( t ) = ![*(*), 0(*)] + ^ (?^l^u-%T,) a -^>U {t ,T,)). (2.32) 
For our harmonic oscillator, the time evolution of the canonical variables is given by 

g(t) = i[if(t),g(t)] = ^|, (2.33) 

p(t) = i [H(t), p(t)) = - p(t) u 2 (t) q(t) , (2.34) 
and by eliminating p(t), we obtain the differential equation 

±( P (t)±q(t))+p(t)u 2 (t) q (t)=0. (2.35) 

This is certainly the equation of motion derived from the Lagrangian L(q, q, t) = p(t) q 2 /2 — 
p{t)u 2 {t)q 2 /2. Note that if we had used W(t,T s ) in Eq. fl2T3T]) instead of U~\t,T s ), the 
equation of motion could not be reproduced correctly when u(t) has an imaginary part. 
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In the Heisenberg picture, the annihilation and creation operators become 



a(t) = U-%T.)a.(t) U(t,T.) = \/^^-q(t) + iJ ^-^-^ pft) , (2.36) 



a(t) = U-\t,T s )a s {t)U{t,T s ) = ^ P -^^q(t) - , (2.37) 
and satisfy the commutation relation 

[a(t),a{t)] = l. (2.38) 
Note that they satisfy the following equations: 

a{t)\0 t ) = 0=(0t\<f{t), ]fr\a(t) = = a\t)]oJ. (2.39) 
The Hamiltonian is then expressed as 

H(t) ee U-\t,T s ) H s (t) U(t,T.) = 0| + P{t)uj2 f q2{t) = u{t) [a(t) a(t) + i] . (2.40) 

We denote by {f(t),g(t)} a pair of linearly independent c-number solutions of (12.35 p . 
One can easily show that their Wronskian, 

W[f,g](t) = f(t)g(t)-f(t)g(t), (2.41) 

satisfies the equation 

^(pW[f,g}) = ff t (pg)-^(pf)g = 0. (2.42) 

Thus the combination (to be called the weighted Wronskian) 

W p [f,g]=p(t)W[f,g}(t) (2.43) 

does not depend on t. Since q(t) is also a solution of Eq. (12.351) . we can expand canonical 
variables q(t) and pit) as 

q(t) = Cl f(t) + c 2 g(t), (2.44) 

p(t) = p(t) q(t) = p(t) [ ci f(t) + c 2 g(t) } , (2.45) 

where C\ and c 2 are some time-independent quantum operators. Equations (12.441) and (12.451) 
can be solved with respect to c\ and c 2 as 
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This can be further rewritten by using (I2.36P and (12.371) as 

' Cl ) = c(t) ( a(tf 

c 2 J \a(t) / 



where 



C(t) 



v(t) v(t) 



W p \f,g]^2p{t)u(t) \-u(t) -u(t) j 

with 

u ® )=p(t)[f(t)±Mt)f(t)], \ v( f\ \= P {t)[m±Mt)g{t)] 

u(t) ) I v[t) 



Note that 



(uv-vu)(t)=2ip(t)cu(t)W p [f,g], 
detC(t) = w ^ ] - ] (= const), 



c-\t) 



-u(t) -v(t) 



y/2p(t)u(t) \ u{t) v{t) 

2.4. Bogoliubov coefficients for finite time intervals 

The Bogoliubov coefficients from time t' to time t are defined by 

'a(t)\fa -p\ fa(t'f 
a{t)) \-P a) ' \a(t') 

Since the operators {ci,c 2 } in (12.47)) do not depend on time, we have 

c 2 J \a(t)J \a(f) 

from which we find 

a ~^)(t-t') = C-\t)C(t') 
—p a ) 

-i (u{t')v{t) -v{t')u{t) u{t')v{t) -v{t')u{f) 

2W p [f,g] ^/pjt)u(t)p(t')u(t') L(t') u{t) - u{t') v(t) v(t') u(t) - u{t') v(t) j 
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This is the fundamental formula to express the Bogoliubov coefficients in terms of a given 
set of independent solutions {/(£), g(t)} . 

Due to the commutation relations (I2.38p . the Bogoliubov coefficients should satisfy the 
relation 

(aa-(3(3)(t;t') = l. (2.56) 
This can be directly checked by using the identity ( I2.5ip as 

(«« - ftp) (t; t') = det C~ l it) det C(t') = = 1 • (2.57) 

** put 91 

Note that a ^ a* and j3 ^ P* when u(t) has a nonvanishing imaginary part. 
2.5. Wave functions 

Using the Bogoliubov coefficients, we can express the Heisenberg operators q(t) with the 
creation and annihilation operators at a different time tj as follows^ 

:(a(t; tj) aj — j3(t; tj) aj + a(t; tj) aj — P(t; tj) aA 



= <p(t; ti) ax + (p{t\ h) a/ , (2.58) 
where we have defined functions tp(t]ti) and <f(t;tj) (to be called wave functions) as 

^ *') = /o L m ("^ *') ~ W> = un f 1/9 — ( Vl f{t) - Ul 9{t) ) ' (2 - 59) 

y/2p(t) u(t) W p [f, g] V2p/ ui 

<p(t; tj) = 1 (a(t; h) - p(t; t T )) = 1 {vj f(t) - u T g(t)) . (2.60) 

y/2p(t)u(t) W p [f, g] V^Pi^i 

By using the t-independence of W p [f, g] and Eq. f 1 2 . 5 [) . we can show that <p(t] tj) and <p(t; tj) 
are normalized as 

W P b(t;tj),^(t;tj)] = i (Vt , Vtj) . (2.61) 
Moreover, the following relation holds for any value of t: 

« -P\ u. t , = (-w p [(p{t;t x )Mt'M -wMWi)Mt-M\ (2 62) 
-/3 a J 1,0 1 1 wXMi),¥>(f;fo)] Wpb(*;*0^(*;*o)] 



8 In the following, we will use the shorthand notation such as // = /(ij) or // = f(tj) when a quantity 
is evaluated at time tj . 

12 



We are now in a position to make a few comments, 
basis independence 

We can show that the Bogoliubov coefficients and the wave functions {(p(t; ij), tp(t; tj)} 
do no depend on the choice of a pair of independent solutions {f(t),g(t)} , as they should. 
In fact, suppose that we take another pair {f'(t),g'(t)}. They should be expressed as linear 
combinations of {f(t),g(t)} of the form 

(f'(t) g'(t)) = (f(t) g(t)) S (He GL(2, Q) , (2.63) 

from which we have 

C\t) = S" 1 C(t) . (2.64) 

The new Bogoliubov coefficients associated with the choice {f'(t),g'(t)} then become 
a' -/T 





P (t;tj) = ic\t)Y 1 c\t I ) = c-\t)~~^c{t 

—p a! 



I 




C- l (t)C(t I )=[ n " )(t;t I ), (2.65) 



which shows the basis-independence of the Bogoliubov coefficients. The wave functions 
{(p(t; ti), <f(t; tj)} are also basis independent since they are expressed by the basis-independent 
Bogoliubov coefficients [see Eqs. (I2.59P and f)2.60p . 

lapse independence 

We can show that the Bogoliubov coefficients and the wave functions ip(t; tj) behave as 
scalar functions under the temporal reparametrizations preserving the foliation of spacetime. 
In fact, for such reparametrization t — > t = t(t) , the pull-back of the lapse function N(t) 
[see Eq. (12. 2p ] is given by N(t) — > N(t) = (dt/dt) N(t(t)) , and we can choose a new pair of 
solutions {f(t),g(t)} as f(t) = /(t(t)) and g{t) = g(t(t)j. Then, we can easily show that 
the functions p(t)u)(t), u(t), v(t), and W p [f(t),g(t)] transform as scalar functions under 
the reparametrization. Since the Bogoliubov coefficients and the wave function ip(t; tj) are 
written as combinations of these functions, they also transform as scalar functions. This 
means that there is no need to care about the temporal reparametrization [i.e. the choice of 
the lapse function N(t)] when we construct vacua. 
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2.6. Feynman propagators 



We consider the region ti < to < {t, < t\ < tf , where tf and ti are the future and the 
past boundaries of the spacetime region we consider. The in-in and in-out propagators are 
defined with the following two steps. 

step 1 

We first introduce the following two-point functions from our wave functions tp(t; tj) and 

£(t;tj)9 



Goo(t, t'\ to, to) 



(0t 




?(*<) 


Oto) 


(o to 


o* > 



W p [<p*(t;t ), <p(t;to)](T a 



£*(t >; t )£(t<;*o), (2.66) 



Gio(t, t ; ti, to) 



(0 h \Tq(t)q(t>)\0 ta ) 



i 



<f{t>;ti) v?(t<; t ) (t: arbitrary) , (2.67) 



where t> = max(t,t') and t< = min(t, t'). 
step 2 

We then define the in-in and in-out propagators by sending to and t\ to the values at the 
temporal boundary: 

G Win (t; ^ = lim Goo ( t t ' ; t0; to) s G in /° ut (t, = lim G 10 (t, ti, to) • (2.68) 

The last expression of (I2.66P is obtained from the identities 

_j W p [<p*(t; to), (p(t; t )}(T s ) f i 



ty p [^(t;t ),^(t;t )](T s ) u W p [^(t;to)^(t;t )}(Z 



(0*0 




a 


0*0 ) 


(0*0 


0*0 


) 



1 



W p [<p*(t-,to)<p(t-,to)]{T a 



(2.70) 



9 If we instead use G' 1Q (t,t';t 1: to) — (0 tl | g^(f>) |0t o )/(0 tl |0 to ), then the corresponding in-out prop- 
agator will not coincide with the propagator obtained by the standard path integral (see section [5]) , and 
thus we do not consider this choice in this paper. By contrast, the in-in propagator still has options for its 
definition (e.g., G' 00 (t, t'; t 0} t ) = (0j | Tq(t) q(t') |0t o )/(0t o |0t o }), and we leave as a future work the detailed 
study of such options as well as the relation to the path integral based on the Schwinger-Keldysh formalism 

mm- 
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which can be shown by using the hermiticity at the time T s , q\T s ) = q(T s ) andp^(T s ) = p(T s ) 
(see Appendix [SJ. We then have 

(°*ol (<p(*>;*o)a + ^(t>;to)ao) (v3(t<;to)a + ^(t<;to)a ) |0 to ) 



Goo(t, t'; to, to) 



¥>*(*>; to) <p(t<;t ) 



(OtoJ Qq Qq |0t ) 

(OtolOto) 



1 



-— -^(t>;t )^<;to). (2.71) 
(*;*o), <p(t;t )](T s ) 

On the other hand, the last expression of (12.671) is obtained by using the following identities 
which are direct consequences of Eqs. (12.531) . (I2.56P and (I2.62[) : 

ai = a(ti, t ) a - (ai + /?(ti; t ) ao) 

(a -/S(ti;t )ai) , (2.72) 



(0 tl | ai a |0 to ) 1 (O tl |(ao-/S(ti;to)ai)ao|0to) 



(OtJOto) «(*i;*o) (0 tl |0 tQ ) afat 



o; 



1 



(t: arbitrary). (2.73) 



We then have 



^ , , ,/ , . v _ (0 tl | (<^(t>; ti) ax + v?(t>; ti) ai) (<p(t<; t ) a + <p(t<; t ) a ) |0 to ) 

^lO^, t ; Cl, t0j — -y— r 



<p{t>;ti) <^(t<; t ) 



ti) ^(t<; t ) (t: arbitrary) . (2.74) 



W p [^(t;ti),^(t;to)] 

When p(t) and u(t) are asymptotically constant in the remote past [i.e. p(t) ~ p m and 

cu(t) ~ Win as t — > ti], we can choose a pair of independent solutions {f(t),g(t)} as those 
which behave as 

f(t)~e-^\ (7(t)~e +i ^ (t~t 4 ). (2.75) 

If we choose such basis, we then have 

u ~0, u ~-2i/> in w ill e- ih «» t0 , u ~ 2ip in w in e iw -*° , v ~ (t ~ ti) , (2.76) 

and from Eqs. ( 12. 59ft and (12.601) the wave functions at the remote past are found to behave 
as 

~ —^e-^ (*-*>>, £(t ; t )~__^ e +i<*» ('-*>) (t ~ t ~ *i) ■ (2-77) 

V ^Pin W in V ^Pin <^in 
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A conclusion of the same kind can be obtained for the wave functions (</?(i;ii), (p(t;tx)) 
if p(t) and u(t) are asymptotically constant at the remote future. This behavior will be 
directly seen in concrete examples given in section [3] and in Appendix [B] 



3. Simple example: scalar field in Minkowski space 

In this section, to demonstrate how the prescription of the previous section works, we con- 
sider a free real scalar field living in Minkowski space with the metric 

ds 2 = -dt 2 + dx 2 . (3.1) 

Another well studied example is investigated within our framework in Appendix [Bl 



3.1. Setup 

In order to clarify the structure of mode functions, we first assume that the spatial part is a 
(d — l)-dimensional torus of radius L/2it , which we will take infinite afterwards. The wave 
vectors k then take the following values: 

k = ^n (neZ"). (3.2) 

For k = (ki, &2, • • • , kd-ij , we write k > (or k < 0) if the first nonvanishing element in the 
sequence {ki, &2, • • ■ } is positive (or negative). Note that k < is equivalent to — k > . 

We introduce a complete set of (real- valued) eigenf unctions {Y kta (x)} of the spatial 
Laplacian A d _x = J2i=i ®1 as 

fc = 0: Y k=0 , a=1 = ^= {V = L d ~ l ), (3.3) 
W 



k>0: Y k>a=1 (x) = cosfc • x, Y k>a=2 {x) = J ^ sin k ■ x . (3.4) 
They satisfy the orthonormal relations, 

/d"*y fc .(*)*.,(*) = 4,«M>, (3-5) 

and we expand a scalar field 4>(x) as 

4>(x) = 4>(t, x) = Y,Yl ^*.«(*) Yk >*W ■ ^ 3 - 6 ) 

fc>0 a 
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The action then becomes 

S[<p{x)] = [ d d x \ - - d"<l>(x) dMx) - r 



2~ 

2 



EE [^[l^-^tijt)], (3.7) 



fc>0 a 



where 



d_1 1/2n 



= Vk 2 + m 2 -ie = oo k -ie (k = \k\ = kf^ ) . (3.8) 

i=i 

We thus have the following correspondence with the ingredients of the previous section: 

q(t) = 4>k,a(t) i pit) = 1 ? w(t) = 0Jk, e (= constant in t) . (3.9) 

3.2. Propagator for each mode 

The equation of motion is given by q + u 2 q = , and we choose a pair of independent 
solutions as 

f(t) = e~ lut , g(t) = e i » t . (3.10) 

Their Wronskian is given by W p [f, g] = W [f,g} = 2 icu , which is constant in t . 
The functions u(t) and v(t) are easily found to be 

{"® }=p[f(t)±iu J f(t)] = { , (3.11) 

I u(t) J I -2iue~ 1U3t 

Up[y(f)±iwy(f)] = ^ n , (3.12) 

I v{t) ) I 

and using (12.591) and (I2.60P we obtain the wave functions as 

tp(t;t ) 



W p V2p uj 
^ = w /o 



h /(*) 


u g(t)] = - 


1 e -iu(t-to) 


(3.13) 


[«b /(*) 


u g(t)] = - 


1 e iw(t-to) 


(3.14) 


[«i /(*) 


uig{t)\ = - 


1 e -ia;(t-ti) 

v/2w 


(3.15) 


[vi fit) 


uig(t)] = - 


1 e iw(t-ti) 

v 7 ^ 


(3.16) 
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The Bogoliubov coefficients then can be calculated by using ( I2.62p as 

a(t 1 ;t ) = e i "^- t °\ P(t 1 ;t ) = 0, a(h;t ) = e^fr-*) , P(t 1 ;t ) = 0, (3.17) 



which indicates that the vacuum at a later time, |0 tl ), coincides with the vacuum at an 
earlier time, |0t ) , up to a phase. The two-point functions take the forms 

(0tol4,a(*>)<^,a(*<)K) _ i£*(t>;to)£(*<;*o) 



Gk,oo(t, t ; to, to) 



(ojo fo ) M^[^(t;to)^(t;to)](r s ) 

e -K, e *> e i^, s *< ; (3.I8) 



2Rew fc>£ e- 2eT < 



kfc, io(t, t ; ti, t ) = == — = . , <£(£>, ti) y?(t<, to J 

(OtJOt,,) a(*i;*o) 

e - iw *. «(*>"*<) , (3.19) 



2u k 

where we have restored the subscript of u, wt E = - ie. Note that the dependence on 
to and ti totally disappear in G kt ij(t,t']tj,tj) for Minkowski space, and thus we need not 
to take the limit to — oo , t\ — > +oo to obtain the in-in and in-out propagators. We find 
that the in-in and in-out propagators coincide after we take the limit e — > , 

G i * /m (t,t f ) = G i » /ou \t,t'), (3.20) 

which we denote in the following discussions by 

G k (t,t') = J_ e -fc» (*>-*<) (=G- k (t,t'j) . (3.21) 

3.3. Propagator in spacetime 

Once the propagators are obtained for each mode, the propagators in spacetime can be 
obtained by summing them over the modes (k, a) . The manipulation is known very well for 
Minkowski space, but we here review it briefly for later reference. 

The in-in or in-out propagator is given by the following summation (ti = — oo, tf — oo): 

<0*ol4a(*>)^V(*<)l°*>> 



G^(x,x') = E 



G in/0Ut (x X')= V V ^llZ^tl^jl ^',a'(t') \0to) 
k,k'>0a,a> \ u *il u W 



Y k>a (x)Y k , ia ,(x'), (3.22) 

to— >ti 

Y k>a (x)Y k ^ a ,(x') . (3.23) 

t^tf 
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Since the two-point functions are diagonalized with respect to the modes: 

= G k (t, t') 5 k k > 5 a:a > , 





<«(*>) 


0fc',a'(*<) 


K> 


(0*0 

{0 tl \T(/> kta (t) 


0*0 ) 


o*„) 


(0*! 


0t ) 



to— >ti 



tO-**! 
tl->tj 



a, a' 



we have 

f in/in 1 
L in/out J 



X, X 



fc>0 



fc>0 



Here, R k (x,x') = J2 a Y kia (x) Y k<a (x') are easily calculated as 

1 2 

R k=0 (x, x') = — , R k>0 (x, x') = — cos k ■ (x - x) 

and thus we have 



G(x, x') = -\ G k=0 {t, t') + 2j2 G k(t, t') cos k-(x- x') 



fc>0 



^ ^ G k (t, t') cos k ■ (x — x) 

k 

d d ~ 1 k 

t-t Gfe(t, t') cos fe • (a; — a;') 



(3.24) 
(3.25) 



') = J2 S ^,a(*) y M(*') = G k (t,t')R k (x,x'). (3.26) 



(3.27) 



(2tt 

where we have taken the limit L — > oo in the last equality. 

The integration (13.281) can be performed easily (see Appendix [C]) , and we obtain 

1 



(3.28) 



G(x,x') 



(2ir)~ \x - x'\ 2 

m (d-2)/2 



poo 

/ dkk^G k (t,t')J^3(k\x-x'\) 
Jo 

_ 2)7I Ka^ (mVa + iO) 



(3.29) 



(2n) d / 2 (a + i0y 

with a = [x — x') 2 . Here, J v {z) is the Bessel function, and K v {z) is the modified Bessel 
function of the second kind. 



4. Scalar field in de Sitter space 
4.1. Geometry and definitions 

We first recall the geometry of de Sitter space and collect the notation and definitions. 
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d- dimensional de Sitter space dS^ has the topology M x S d 1 and is defined as a hyper- 
boloid 

VMN X M X N = e (M,N, - ■ - = 0, . . . ,d) , 

(vmn) = diag(-l, 1, . . . , 1) , (4.1) 

in (d + l)-dimensional Minkowski space with the metric 

ds 2 = t]mn dX M dX N . (4.2) 

£ is called the de Sitter radius. The constant Ricci scalar curvature is then given by R = 
d{d-l)/£ 2 . 

There are several well-known coordinate patches which cover the whole or just a part 
of de Sitter space. Among them, we consider only the global patch and the Poincare (or 
planer) patch, which we will briefly review below. 

global patch: This coordinate patch covers the whole region of de Sitter space. The 
embedding of dS^ is given by the following functions: 

X°(t,Q) = £ sinhr, X 7 (r, Q) = £ cosh r tt 1 (J = l,...,d) (4.3) 

with f2 • O = 1 . Here, r runs over the range — oo < r < oo and f2 is a unit vector in M. d 
spanning a (d— l)-dimensional sphere. With the coordinates (r, ft) the metric has the form 

ds 2 = ^ [ _ dr 2 + cosh 2 rd ^2_J = f[ _ ^ _ t 2y2 dt 2 + t 2yl dQ 2_J (44) 

In the last equality, we have introduced another temporal coordinate t as 

t = tanhr (-1 < t < 1) . (4.5) 

Poincare patch: This coordinate patch covers only half of de Sitter space. The embedding 
is given by the following functions with rj < and x G 



pd-l. 



X°( V ,x) = e V \ +IX] \ X\ V ,x) = £^, X d ( V ,x) = f + v2 lxl \ (4.6) 
—It] —rj —2r] 



where the spatial norm is defined by \x\ = y/x ■ x = y/ 5ij x % xi . Note that this patch only 
covers the region X° + X d = £ 2 / {—rj) > . In these coordinates, the metric takes the form 

d, 2 = £ 2 ~ ^ \ dx ■ dx . (4.7) 
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We define the de Sitter invariant quantity 

Z(x, x') = r a vain X M {x) X N (x') , (4.8) 
which is related to the geodesic distance d(x, x') between two points x and x' via the relation 

Z(x,x') = cos(^^-). (4.9) 



It takes the following values depending on the positional relation between x and x' : 

> 1 (for x and x' timelikely separated) 
Z{x,x') \ =1 (for x and x' lightlikely separated) , (4-10) 
< 1 (for x and x' spacelikely separated) 

as can be seen from the identity 

Z(x,x') = 1 - ^ mN {X M {x) - X M {x')) (X N (x) - X N {x')) . (4.11) 
In the global and the Poincare coordinates, Z(x,x') is written in the following form: 



- sinh r sinh r' + cosh r cosh r' Q ■ Q' (global coords.) 
ZM = < v 2 + r] '2_ lx _ x f (4.12) 

(Poincare coords.) 

27]7]' 

Any de Sitter invariant quantity depending only on the geodesic distance d(x, x') (such 
as Green's function of a de Sitter invariant vacuum) should be a function of Z(x,x') , and 
actually all the propagators considered in this paper can be written as functions of Z(x, x') . 
In what follows (except in subsec. I6.3p . we set t = 1 . 



4.2. Scalar field in the Poincare patch 

We first consider a free real scalar field in the Poincare patch. The action takes the forrrj^ 
S e [<Kx)] = - \ J d vJ d^x^ig^d^d^+im 2 - ie) 2 ) . (4.13) 

Using the same eigenfunctions {Y k a (x)} as those given in section [3] [Eqs. (13. 3p . f l3.4p ], we 
expand a scalar field 4>(x) as 

<l>(x) = <j>( v , ^) = zZzZ foM Y >*A X ) ■ ( 4 - 14 ) 

fc>0 a 



10 In this paper, we put a possible curvature-coupling term, (£/2) R<p 2 = (d(d — 1)^/2) </> 2 , into the mass 
term, (m 2 /2) (j? . 
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The functions defined in section [2] then take the following form (see (12. 8p ): 

q{v) = Mv), P(V) = (-v)- {d - 2) , 

u)(rj) = \J (m 2 — ie) (— r])~ 2 + k 2 (k = \k\) , 

from which we introduci"] 



u = u(r] 



r/o— oo 



-Vi 



mi 



m — ie (for m > 0) 

— rji k — ie (for m — 0) 



(4.15) 
(4.16) 

(4.17) 
(4.18) 



4.2.1. Propagators for each mode in the Poincare patch 



The equation of motion (I2.35P takes the form 

rf q{rj) — (d — 2) q{rj) + (k 2 rf + m 2 — ie) q{rf) = , 
and we choose a set of independent solutions as 

f(rj) = (-r/)^ 1 Jv a (-krj) , g(rj) = (-r/)^ N Ve (-krj) 

with 



(4.19) 



(4.20) 



yj \^) 2 - (m 2 - ie) = v + ie (m < ±£) 



i Jm 2 - ie - (^=i) = i// + £ (m > ^) 



z/ = 



d- 1\ 2 



w? — 



d- 1\2 



(4.21) 



2 / V V 2 

Here, N u (x) is the Neumann function. Note that Re v e > for any positive value of m . The 
Wronskian is given by W[f,g](rj) = —(2/ir) (—r]) d ~ 2 , and thus 

W p \f,g]=p(t)W\f,g](t)= --. 

71 

The functions u(rj) and v(rj) of (I2.49P have the form 



(4.22) 



u{r}) 
v(rj) 



-rj) 

-rj) 



d-l 
2 



d-l 

2 



d- 1 
2 

d-l 



i/ E ± itu 77 J r]) + kr] J x+Ue (-k 77) 

u e ±iur}) N Ve {-krj) + k rj N 1+Ve (-k rj) 



(4.23) 



(4.24) 



11 Our prescription of section [2] cannot be applied directly to the exactly massless case, where uj(rj) =0 for 
k = . Wc actually define the massless theory as the m — > limit of a massive theory, and the expression 
for fh\ in the massless case is presented just for reference. 
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where we have used the formulae 

z = v J v {z) - z J v+1 {z) , z = vN u {z)-z N v+1 (z) . (4.25) 

oz oz 

The wave functions are then given by 

pfaifr) = --^= }~ V) ]_ 2 (vjJ^-k^-ujN^-kr))) , (4.26) 

Tp( v;Vl ) = (vjJ^i-kri-ujN^-kr))) , (4.27) 
with which we obtain an exact expression of the propagators for a finite time interval: 
G k ,oo(v,v';Vo,Vo) 

-i7r[(-77)(-7/)]^ 



2(u*v - v% n ) 

x[voJv-A-kr/)-uZN v _ s (-krfj\ [v J Ve (—k rj') - u N Ve (—k r]')] (m < — ) 
-ivr [(-//) (-?/)] V 



2[(uq?Jo - ^o u o)cos(7rz/„ e ) + (|n | + |v | ) sin(7rz/_ e )] 

(4.28) 

= Z^Kz!?Kz!^Z: J, E (-fc77) - uiJV* (-fcifl] [tJb J^(-*i?0 - uoJV^C-Ai/)] • 
2(«i w - «i Mo) 

(4.29) 

We now send 770 , 771 to the boundary of the Poincare patch; 770 — > — 00 and 771 — ^ . The 
wave functions 99(77; 777) and (^(77; 77/) (/ = 0, 1) then behave as 

<P(v,Vo) -^M^^) (-rj)^H£\-kri), (4.30) 



<p(mio) ~ ^ e -^ + ^) (-rj^Hlfi-kr,), (4.31) 

<p{m rn) ~ - V{V i^' V£ (-Vi)^ -S- imi) (-77)^ JuA-kri) , (4.32) 

r(u e )(k/2)-"' , „ /d-l \ . .d-i . , . , . 

^(77; 77O 1 (-771)^ (— - ve + imi) (-77)— • (4.33) 

Here, Hv(x) and H^\x) are the Hankel functions defined by H^ ,2 \x) = J u (x) ± iN u (x) . 
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They can be obtained by using the asymptotic forms 

j> ~ ±i\l— (-Vo)-— e ±1 l**»+ — a — J , (4.34) 

(-^-Ve^h+^i, (4.35) 



«i 1 (k/2) Ue . . d -i ,„ \ 

which can be easily derived by using the asymptotic forms of the Bessel functions: 

J *<"*i) ~ TvT^j ( 2 ) ' ^(-^)'~~(-2) ' (438) 

Jj-fcT?) ^-{-k^Y 1 / 2 co S [-k V - (7r/4)(l + 2i/ e )] , (4.39) 

N„,(-kv) n ~^° ^(- k v)~ 1/2 sm[-k7]- (7r/4)(l + 2i/ e )] . (4.40) 
Note that the wave function </?(t7; 770) converges to a finite function, 

<Pm(r,) = ^(-v)^H^(-k V ), (4.41) 

in the limit ?7o — >■ — oo (up to an oscillatory phase), while (p(rj;rji) is divergent in the limit 
?7i — ^ 0. This difference can be attributed to the fact that the timelike vector becomes 
asymptotically a Killing vector in the remote past, but does not in the remote futureo The 
finite asymptotic function (p- m {rf) coincides with the positive- mode wave function associated 
with the Euclidean vacuum up to a phase. 

The Bogoliubov coefficients have the asymptotic forms 

aiVl]r]0) } ^L(-fcW2)--f— -^-i^e^C^^), (4.42) 

where we have used the approximation (— r/i) _Ve ^> (— 77i) 1/s , since z/ £ always has a positive 
real part. Note that the norm of the Bogoliubov coefficients 1 01(771; f]o) | and | /3 (771 ; 770) | diverges 
in the limit 771 — >■ although the relation ( 12. 56ft still holds. This is an evidence that the 
remote future cannot be treated adiabatically [TJ. 



12 In fact, the Lie derivative of with respect to the vector £ = d v is £^ g^ u oc (—rj) 3 
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Substituting everything into Eqs. (I2.66P and (I2.67P and taking the limit r] — > — oo and 
r]i — > , we obtain the propagators in the Poincare patch, 

Gt /m (v, V) = \ [(-vX-rf)]*? HU(-k V> ) H%(-k V< ) , (4.43) 
Gt /0U \V,V') = \ [(-vX-rf)]*? J^-kr)>) H^(-k V< ) . (AAA) 

One should note that, in the limit 771 — > , the divergence in ip(r]]r]i) is canceled out with 
that of a (771; 770) and the in-out propagator is obtained with a finite value. 

4.2.2. Propagators in the Poincare patch 

Since the eigenfunctions {Yfc ja (a:)} are the same as those given in section [3], the propagators 
in spacetime can be written in the form 

fc>0 a 



d-1 , 1 d-3 

f/\ 2 



(2ir)^~ \x — x 
as in the case of Minkowski space (see Appendix |C]) . 
For the in-in propagator, we have 

G' m/in {x,x') 



POO r in/in i 

/ dkk^ J^(k\x-x'\)G k in/out (r),r)'), (4.45) 
Jo 2 



^ (27T) 2 \x — x' 



roo 

\ dkk^ J i ^(k\x-x'\)Hj > 1 ] e (-kr] > )Hlf(-krj < ). (4.46) 
Jo 2 



As is proved in Appendix |Dl this can be integrated to the form_f 



2 (27r) d / 2 



d-2 



with u = —Z(x, x') + i . This can be further rewritten as 

GWfc z<) = + F (<L± + v ,izl- v .*.l-±\ (4.48) 

V ' ; (4vr) d / 2 r(d/2) V 2 ' 2 ' 2 ' 2 / V ; 

rc^=±i d-i 

(TV, («) , (4.49) 



4tt(^ 1 )/ 2 sin[7r(^fi-r/)] 



13 Here we have set e = . V£(z) and Q^(z) denote the associated Legendre functions of the first and 
second kind that are defined on the complex z-plane other than the cut along the real axis to the left of the 
point z = 1 . There are another type of associated Legendre functions that are defined on the interval (—1, 1), 
which we denote by P(a;) and Q(x) . See Appendix lEl for their definitions and several useful identities. 
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where Rummer's relation (1B.29|) has been used in the first equality and C£(x) is the Gegen- 
bauer function. This propagator is the same as the well-known in-in propagator associated 
with the Euclidean vacuum. 

The in-out propagator 
G' m/out {x,x') 



2 /~ \ d ~ l i .1 d ~ 3 



(2ir)— \x - x'\ 2 Jo 
can be rewritten in a similar manner to the form 

-i7r(d-2) 



poo 

/ dkk^ J^B(k\x-x'\)j Ue (-kri > )Hlf(-kT ]< ) (4.50) 
Jo 2 



G^\x,x') = (u 2 - 1)-^ QZi /2 (u) (4.51) 



with u = Z(x, x') — iO . A proof is given also in Appendix |D] 

If we consider the massless limit where v — > (d — l)/2, the in-in propagator diverges, as 
was pointed out in [3]. In contrast, we find that the in-out propagator has a finite massless 
limit: 

G in /° ut (x, x') = e ( J^ /2 2) (u 2 - l)-^ 2 sS(n) (massless) . (4.52) 



4.3. Scalar field in the global patch 

In the global patch, we choose a set of (real- valued) spherical harmonics {Yij^(ft)} as the 
complete eigenfunctions of the spatial Laplacian A^-i on S^" 1 . They satisfy A^-i Ylm = 
—L{L + d — 2) Ylm (M = 1, . . . , N^) , and the degeneracy is given by 

N ^ = % + -2)m {2L + d - 2) (4 ' 53) 

with the exception of the case d = 2 and L = where Nq = 1. We choose them such that 
they are orthonormal: 

y dft y LJlf (n) *W(n) = $ll> Smm' ■ (4.54) 
Then, by expanding 4>{x) as 

0(x) = ^0 iM W^LM(^), (4.55) 

L.M 
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the mode function q(t) = 0la/(^) describes a harmonic oscillator with time-dependent mass 
and frequency of the following form [see (12.81) ]: 

p(t) = (l-t 2 )^, (4.56) 



(t) = (l-t 2 ) 1 Jl (L + d - 2) (1 - t 2 ) + m 2 - is = (l - t 2 ) 1 fh(t) . 



(4.57) 



4.3.1. Propagators for each mode in the global patch 



We choose a pair of independent solutions as 



/(*) = (1 " t 2 )^- P?{t) , git) = (1 - t 2 )^ Q*(t) 



where 



k = L 



d-3 



(4.58) 



(4.59) 



and P£(i) and Q^(£) are the associated Legendre functions defined on the interval (—1,1) 
(see Appendix [Ej) . The weighted Wronskian (constant in t) then has the form 

r(k + u e + i) 



W p [f,g]=p(t)W[f,g](t) 



T(k -u e + l) 

The functions u(t) and v(t) defined in ( 12. 49ft are then given by 

u(t) 
u(t) 

v(t) 
v(t) 

and the wave functions (if, (p) take the form 

r(k -v e + i 



(4.60) 



(l-t 2 y * [(Lt±im(t))P£(t)-{k-v £ + l)PZ +1 (t) 
(1 - t 2 )^ \(Lt± ifh(t))Q£(t) - (k - v e + l)Q^ +1 (t) 



<p{t\ti) 



(l-t 2 )^( Vl f{t)-u ig {t)) 



y/2mj T(k + v e 

T(k — v E + 1) . 9 . d-i , „. . , . N 

/~— r(h l Lu ^ - # ~ fa / * - • 

y/2m I T{k + i/ £ + 1) 



(4.61) 
(4.62) 

(4.63) 
(4.64) 



We now send to ; ti to the boundary of the global patch; t$ —1 and t\ — > 1 . Using 
Eqs. (IE.21j) and (IE.22j) . we see that iti = u(ti) and t>i = v(ti) take the following asymptotic 
forms in the limit t± — > 1 : 

d- 1 



"i 

Ml 



2" e sin(vrz/ £ ) r(v e 



7T 



~ 2^- x cos(7rz/ e ) r(i/ e ) 



2 

d-1 



d— 1 U£ 



+ v £ ±im 1 ) (1-tj) 4 



(4.65) 
(4.66) 
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and thus we find that the wave functions behave as 14 



tp(t\ti 



2^- 1 T(u P )/ d-1 



y/2mi 



+ u £ ± imA (i - 1\) 



2 



(l-t 2 )^P"^(t). (4.67) 



Here, we have used the fact that mi ~ m — is in this region. 

Similarly, no = u(to) and vq = v (to) have the following asymptotic forms in the limit 
to -> -1: 

( 2~ v ' cos(irk) T(k + v £ + l) { d-1 

u } 

u J I 2"- sin(rt) T(^) f d - 1 _ , _ , 21 -i=i-f 



r(„ e + i)r ( *-, e + i) v 2 "^Nt'-'Sr" ( d:odd » 
^«K-r*)rC.)^-i ,, ±im „) (l _ ia 



7T 



i'f) 
^0 



2" e_1 cos(vrA;)r(z/ £ 



2 

d-1 



- z/ e ± im (1 - t 2 ) 4 



d-1 i^e 
2 



2-^- 1 7 rr(fc + v £ + 1) /d-1 
sin(7rA;) T(z/ £ + 1) r(fc - v £ + 1) V 2 



±im ) (1-to)' 



4 ' 2 



(d : even) 
(4.68) 

(d : odd) 

(d : even) 

(4.69) 



Here, we have used the relation (l — to) 2 3> (l — to) 2 > which holds since Rev £ > 0. We 
thus find that the wave functions behave as 



<p(t;t c 
<p(t;t c 



2 V *~ X cos(ttA;) V{v £ ) T(k - v E + 1) ^d - 1 
v^r^r(A; + z/ £ + l) 



±im ) (1 — tg)" 



x 



(l-^)^P^(t) 



(4.70) 



in odd dimensions, and 



<p(t;t ) 

<p(t;t ) 



2 V * sin(7r£;) T(u e ) T(k -u e + l) (d-1 
■K^/m^V{k + u e + l) 



±im ) 



x 



(1-t 2 )" 1 Q?(t) 



(4.71) 



in even dimensions. 

From (I2.62p . the Bogoliubov coefficient a(ti;to) can be found to have the asymptotic 



14 The asymptotic forms given in Eq. (|4.67[) do not satisfy the normalization condition (|2.61l) . We actually 
need terms proportional to (l — t^)"^ 2 , which are omitted in deriving the asymptotic forms, in order to 
ensure the condition. However, these asymptotic forms are yet sufficient for calculating various propagators. 
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form 

a(ti,t ) 



2i sin(7r^ £ 

71 
w E -l 



2^- 1 T(v £ )f d-1 



+ u £ + ifh 1 ) (l-tj) 



x 



cos(Trfc) r(i/ e ) T(k - v £ + 1) /d-1 



v/MoT(fc + i/ e + 1) 



z/ £ - im ) (l - t 2 ) 



(4.72) 



/5(ti;t 



2i sin(7rv F 



7T 



jit— i 



I> £ )/ d-1 



x 



2^~ 1 cos(ttA;) r(z/ e ) r(ife - v e + 1) /d - 1 



^MT Q T(k + v £ + 1) 



im ) (l - *q)" 



(4.73) 



in odd dimensions, and 



Oi(tl]to) ~ — i cos(7ri/ £ ) 



ii/ e -i 



r(z/ 6 )/ d-i 



\J2m\ V 2 
2" e sin(vrfc) T(u e ) F(k -u £ + l) /d-1 



+ i/ B + imi ) (1 - t\) 



2 



7T 



V2W^r(fc + z/ £ + i) 



u £ - im ) (l - ig) 



He 



(4.74) 



/3(*i;*o) 



2i sin(7rz/ £ 



7T 



)i/ £ -l 



r(i/ e )/ d-i 



A/2m7 



imi)(l-^) 



x 



2"» sin(Trfc) T(z/ £ ) r(Jfe - i/ e + 1) /d - 1 



imoj(l-tg) f 



(4.75) 



7rV2^or(fc + v e + 1) 
in even dimensions. 

Although cp(t; ti) and (p(t; t ) asymptotically diverge as (l — tf) 2 and (l — £q) 2 , 
respectively, in the limit t\ — >■ +1 and to ~~ ^ — 1 j we still obtain a finite result for the in- 
out propagator since the divergence cancels out with that of a(ti, to) (as in the case of the 
Poincare patch): 

1 



G^ /out (t,t') 



lim 

*o- ; i a(t i; t ) 



ITT 



2 sin tiv f 



[(i - ti) (1 - 1 2 <)] (d - 1)/4 pr E (^>) p^(^<) (d ■■ odd) 



— — [(1 - ti) (1 - tl)f- 1)/4 P^ (t> ) Q*( f< ) (d : even) . 

COS 7Tl^£ 

Once this expression is obtained, we may take the limit e — > 0. 
On the other hand, the in-in propagator takes the form 

G^(t,t') 



(4.76) 



lim 



to^-iW p [<?*(t; t ), <p(t;t )](T a ) 



t ) *o) 
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[(1 - tl) (1 - t\)] ^ Pf{ t> ) P?(f<) (d : odd) 



- Ts ^ [(1 - tl) (1 - tl)] Q? (t>) Q£(t<) (d : even) . 

I W p [Q k s (t) : Ql e (t)](T s ) >n <Ji hK>) kK<) K > 

When the mass is large (m > (d — l)/2; z/ e = iyU + e), the in-in propagator has a finite limit 
as e — > : 

7T 



[(1 - * 2 >) (1 - t 2 <)] ^ PT(t<) (d : odd) 



G^ m (t,t')=l 2sinh 2 M L ' (4.78) 

1 ^m>) t (1 " ^ (1 " ^ ~" Q ^ (t>) Q " (t<) (d : even) • 

In contrast, when the mass is small (m < (d — l)/2; v e = v + ie), the Wronskian in the 
denominator of (I4.77P vanishes as e — > 0, and the in-in propagator diverges. 

The wave functions at the remote past and future had been obtained for the heavy mass 
case (m>(d- l)/2) in $\5\ afl 

^ out (t)oc(l-i 2 )^P^(t), (4.79) 

VinW oc <^ „_! 4.80 

I (1-t 2 )— Q^(t) (cheven) 

by requiring that <fi* n (t) (<f ut(t)) is regular for t — > — 1 (t — > +1) and an analytic function 
in the lower half of complex m 2 plane (see also |13j where they are obtained by suitably 
choosing the Jost functions). Our propagators (14.761) and (I4.77P for m > (d — l)/2 are 
consistent with these wave functions. 



4.3.2. Propagators in the global patch 

We now take the summation over all modes to obtain the de Sitter invariant propagators, 
G m ' ont (x,x') and G in//m (x, x'). For d > 3, the summation over M can be written with the 



15 We here give the following identities which are useful in comparing our results with the literature: 
2^ +L cosh L Te^^- L ^ T /d- 1 d-1 



r(i + u) 



F ( —Tr 1 + L i —?T~ + L + u;l + v; -e~ 2 A = (l - t 2 ) P7 v (t) ; 



2T i + i cosh L re(V+ i -^)^ /d-l d-1 n 2 ^ ™, n 



r(i - v) V 2 '2 



r(fc + i + ^) e j (-i) fe Pr(*) ( d:odd ) 

r(fc + l-^) 1 (-lJfe+Va (2/^-)Q-"(t) (d:even) 
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Gegenbauer polynomials: 



N 



(d) 



g Y LM (Q)Y LM (Q') = {d _ 2) ^- ^ C L 2 (n.flO (|n i _ 1 |=27r»/r(d/2)). (4.81) 
As for d = 2, the sum has the form 

f - (L = 0), 

E^)«= 2j ■ (4.82) 

m=i — (L>1), 

which is the same as the d — > 2 limit of the expression (14.81]) . Thus Eq. (I4.8ip can be 
understood to hold for any dimensionality d > 2. The in-out propagator in spacetime then 
takes the form 

L=0 M=l 



:0 (d-2)|nd 



which becomes 



2 sm(7rz/) (a — 2) |ii<j-i| 



oo 

X 

L=0 



£(2L + d - 2) p-^>) P£(t<) C L " 2 (cos£) (4.84) 



in odd dimensions, and 

cos(7rz/) (a — 2) \ild-i\ 



oo 

X 

L=0 



£(2L + d - 2) p-"(t>) Q£(t<) C L " 2 (cos£) (4.85) 



in even dimensions. Here, we have defined 9 via the relation f2 • f2' = cos#. 
Using Eqs. ( IE. 251) and ( IE. 261) and introducing 

u±(x, x') = - Z(x, x') ± iO = ~tt +cos9 ± iQ 

(1 - t 2 )* (1 - t' 2 )* 
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we can rewrite the in-out propagator in a de Sitter invariant form: 

(ul - lY^r Q7; („ ) - (n 2 „ - 1)-^ Q% U 



^-(in/out/ 
^odd \ X i X 



2 (2tt) 2 sin(7r^) 



le _ 17r _ 



(4.87) 



2 (27r) 2 COs(7TZ/) 



« - 1)"^ s;.K) + (ui - 1)-t- 

(4.88) 

In the massless limit m -)■ (or f — )■ (d — l)/2), we have 

sin7r^ — 7-0 (d : odd) , cos7ri> — >• (d : even) , (4.89) 

and thus the propagators (I4.87P and (14.881) diverge. We thus conclude that there exists no 
finite massless limit in the global patch, as opposed to the case of the Poincare patch. 

On the other hand, the in-in propagator in the heavy mass case (m > (d — l)/2) takes 
the form 

G^(x, x>) = £ f + d ~ 2 Gf in (t, f) Cf (O • IT) , (4.90) 



which becomes 

^in/in 



G - <^'> - 2M^m^< [{1 - tl)<l - tl)] " 



-i 



oo 

X 

L=0 



£(2L + d - 2) p-"(t>) P£(t<) C L " 2 (cos£) (4.91) 



in odd dimensions, and 



ao = . , * 9 . , i [(i - * 2 >) (i - m ^ 

7rsm(7rz/) (a — 2) 



X 

L=0 



£(2L + d - 2) Q-"(t>) Q£(t<) C? 2 (cos 9) (4.92) 



in even dimensions. Here, v — i/i — i-^/m 2 — (d — l) 2 /4 (/i G K+). Note that G™[™(x, x') = 
^odd Ut ( x ' x ')i w hich is consistent with a well-known fact that the in- vacuum equals the out- 
vacuum up to a phase in odd dimensions (see [5]). The summations in ( 14.9ip and ( 14.92p can 
be carried out analytically by using ( 1E.25I) and (IE. 281) . and are again expressed in the de 
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Sitter invariant forms, 



^odd ( x ) x ) — ^odd ( X ' X ) 



in/out , 



;„ d-2 

- e" 17r — 



2 (27r) 2 sin(7rz/) 

^cvcn { x : x ) 



;_ rf-2 
6 2 



2 (27r) 2 sin 7r^ 



>2 -1)-— QjAu, 



[ul-iy— q„ 2 Au, 



+ 



1 7T 



COS 7T^ 



V - 1)"— Q \{U. 



2 (27r) 2 sin 7ri/ 

7T 



[e" 1 - « - 1)-^ V v 2 i + (u 2 „ - 1)"^ P * [ (u. 

"2 v 2 

d-2 _ d ~ 2 „ d-2 _ d ~ 2 , 



COS 7^ 



1)-— p 2 1 (-n_ 



(4.93) 



(4.94) 



5. Feynman path integral in de Sitter space 

In this section, we consider the Feynman propagator obtained by the path-integral, 

j[d<j)}<f){x)<j){x')e iS ^ 



(0(x)0(x')> 
with S e [<f>] 



l - Jd d x^/=g~(g^d^d u <f>+(m 2 -ie)(f> 2 ). (5.1) 



We assume that the base space where 4>(x) lives has a sufficiently large noncompact region in 
the temporal direction near the future and past boundaries at t = tf and t = U, respectively 
(see Fig. [2]). Then, due to the existence of the ie term, when sending the initial time 



t = t 




sufficiently long 
for the projection 
to the ground states 
to be completed. 



t = U 



Figure 2: The spacetime region where the path integral is performed. 



to and the final time ti (in the Hamiltonian formalism) to the infinite past at ti and the 
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infinite future at if, respectively, the initial and final state would be well kept subject to the 
projection to the instantaneous ground state at each moment to or t\ . It is thus natural to 
expect that the propagator defined in (15.11) coincides with the in-out propagator constructed 
in the preceding sections. In the following, we give an evidence for this assertion by showing 
that, for each mode, the propagators numerically calculated from ( 15. ip coincide with our 
in-out propagators^ 

We first expand <fi(x) as in (12.51) : 

<P(x) = Y,Mt)Y n (x). (5.2) 

n 

The propagator is then written as the sum of the propagators of the mode n : 

{<j>{x) <P(x')) = J2(Mt) Mt')) Y n (x) Y n (x') , (5.3) 

n 



f[d(j) n ] e iS ™. s ^»] 



With S n , £ [(f) n } = I di ^ <j) n (t) [-dt opod t - ^(i)] n (i) . (5.4) 



We redefine the variable 4> n (t) by Xn{t) = P 1 ^ 2 (i)0n(i) = 0„(i) and omit the mode 
label n for simplicity. Then, the action for each mode has the form 

Se[x] = Ai \ x(t) [-d 2 t - n 2 E (t)] X (t) , (5.5) 

J ti 

n 2 £ (t) = u 2 £ (t) - (d(i)) 2 - a(t) , (5.6) 
and the propagator for each mode is given by 

(0(i)0(i')) = [p(t)p(t')]- 1/2 (x(t)x(t')} = [ P (t) p(t')]~ 1/2 (t\ _^_ Q2 \t'). (5.7) 

We numerically evaluate the propagator ( 15. 7p by dividing the interval (ij, i/) into N parts 
and by calculating the inverse of the matrix corresponding to i _1 (— <9 t 2 — Q 2 ). We take a 
uniform spacing a = (tf — t{)/N for brevity and write the time variable as i = ar with r an 
integer in the region < r < r ^ (rj = ij/ 'a and rj = i//a). We then introduce dimensionless 
variables Xr as 

X (i) = a 1 / 2 Xr, (5.8) 



16 Note that one does not need to specify the foliation in defining the path integral (j5.1j) . and the obtained 
propagator must be de Sitter invariant as long as the symmetry is not broken by any of the UV and IR 
regulators (which are absent in free field theory considered here). The de Sitter invariance of the in-out 
propagators obtained in the previous section is an expected property if this equivalence is correct. 
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with which the action becomes 

S e[x] = \ ^S^XrXs 



Since 



S e ,rs = [2 — a 2 Q e (ar)j S r>s — S rtS+ i — 5 r>s _i 



, . Jd N x XrXr'exp(i J2S £:SS >XsXs>) ., c _ ly 



Jd N x exp(i 52S EtSsl x s Xs>) 
the propagator is obtained as 

(0(t)0(t')) = ia[p(t)p(t')]" 1/2 K 1 ) rr , (* = ar, t' = ar'). 



(5.9) 
(5.10) 

(5.11) 



(5.12) 



We numerically calculate the inverse matrix f)5.12p for both of the Poincare and global 
patches and compare the result with our in-out propagators obtained in section HJ 

The result in the Poincare case is depicted in Fig. EJ while that in the global patch of 
d = 3 (d = 4) is in Fig. H] (Fig. E]). We find that there is a perfect agreement for the 
global patch, while there exists a small discrepancy for the Poincare patch. We observe that 
the discrepancy gets reduced as one takes a finer mesh near rj = and a larger value for 
1 770 1 , and expect that it will disappear eventually. We thus are almost convinced that our 
in-out propagators coincide with those obtained by the Feynman path integral with the is 
prescription. 





Figure 3: The real part (left) and the imaginary part (right) of the in-out propa- 
gator G)^ out (?7, 7]') (solid curve) and the propagator (4>k{v) 4>k{v')) (dashed curve) 
for d = 3, m = 4, k = 1, a = 0.002, e = 0.8, rj = -20, 771 = -0.01, and 
77' = —10.005. Recall that 77^ = —00 and 77/ = . 
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Figure 4: The real part (left) and the imaginary part (right) of the in-out propa- 
gator G^ out (t, 0) (solid curve) and the propagator (0l(£) <Al(0)) (dashed curve) 
for d = 3, m = 0.5, L = 1, a = 0.005, e = lO" 10 , t = -0.995, and t x = 0.995. 
Recall that t{ = —1 and tf = +1 . 




Figure 5: The real part (left) and the imaginary part (right) of the in-out propa- 
gator G^ out (t, 0) (solid curve) and the propagator 0l(O)) (dashed curve) 
for d = 4, m = 0.5, L = 1, a = 0.005, £ = 10~ 10 , f = -0.995, and i x = 0.995. 
Recall that tj = —1 and t/ = +1 . 
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6. Heat kernel representation and the composition prin- 
ciple 



6.1. General theory 

We consider the random walk of a relativistic particle moving in a Lorentzian manifold with 
the metric 

ds 2 = g^(x)dx^dx u . (6.1) 

Its trajectory is uniquely specified by the functions X fl (X) (0 < A < 1) , up to reparametriza- 
tions A — > /(A) such that d/(A)/dA > 0, f(0) = 0, f(l) = 1 . The amplitude connecting 
two points x and x' is then given by the Feynman path integral: 

A(x, x') = [ X(1)=X } d ^ X J\ e , (6.2) 



^(o^VoKDiffi) 

where Vo^Diffx) is the gauge volume of one-dimensional diffeomorphisms 

X"{\) X"(A) = X»(f(\)) , (6.3) 
and we propose to set the action Io[X] for the random walk in a Lorentzian manifold as 

Io[X(X)] = -(m-ie) / dA J- g(a/ (X(X)) X»(\) X"(X) - ie' . (6.4) 

Jo 

Note the presence of two infinitesimal imaginary parts, ie and ie' , in J [X(A)] (e, e' > 0). 
The first (ie) is the standard one, which manifestly suppresses the contribution from such 
paths that are prolonged in the timelike direction. We further have introduced the second 
one (ie') in order to define the path integral for any shape of path in a Lorentzian manifold. 
In fact, for a timelike segment (X 2 < 0) we can neglect e' and the action becomes the 
standard action for a timelike path, while for a spacelike segment (X 2 > 0) we can rewrite 
the square root as 



V-X 2 - ie' = Ve-^-VX 2 = - iV^ 2 (6.5) 

and the action gives the path-integral weight which suppresses the contribution from such 
paths that are stretched largely in the spacelike direction. 

The action (16. 4p of Nambu-Goto type is equivalent to the following action of Polyakov 
type: 

J[X(A),e(A)] = y o dA [ — g 6 ] ' (6 - 6) 
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where e(A) > is the einbein defined on the one-dimensional manifold. We can see from this 
expression that e and e' give imaginary parts of the same sign. The new action I[X(\), e(A)] 
also has the invariance under the one- dimensional diffeomorphisms 



X"(\) -+ X"(\) = X"(f(\)) , e(A) -> e(A) 



d/(A) 
dA 



e(/(A)) 



(6.7) 



and we can take a gauge fixing where e(A) = constant = T . However, as is discussed in 
detail in [13], such constant T = J^dA e(A) is actually Difivinvariant and needs to be further 
integrated, so that we obtainF^l 



A(x,x') 



X(l)=x 



X(0)=x> 



[dX»(\)de(\)] ihmx)] 
Vol(Diffi) 



dT 



X(l)=x 



[dX(X)} exp 



X(0)=x> 







dTe - £ 7(2T) 



X(T)=x 



[dX(t)] exp 



2T 

T 

dt 



X 2 (t) + ie' m 2 - ie 



(6. 



'X(0)=x' 

where in the last line we have rewritten the expression with t = T A . The path integral is 
nothing but that for the quantum mechanical amplitude from the state \x') to the state \x) 
with the Hamiltonian 

1 



H = - (-□ + m 2 - ie) = - 
and thus we obtain the expression 



-9 



(6.9) 



A{x,x')= / dTe-^'W ( x \ e -Wm-v+m 2 --)| x ') 
Jo 

poo 

= / dT e- 1 ™ (- 2 -i £ )-e 7 (2T) K ^ x ,. T) 
Jo 

Here, K(x, x'\ T) is the heat kernel of the d'Alembertian □, 

K(x,x'; T) = (x\e i( - T ^ D \x') , 
which satisfies the following equations: 

i A K(x, x'- T) = ~D X K(x, x'- T) , 



K(x,x'; T = 0) = 5 d (x,x') 



-.9 



5 d (x - x') . 



(6.10) 

(6.11) 

(6.12) 
(6.13) 



17 There may arise a divergence when calculating the Jacobian to obtain the second line, but such diver- 
gence should be ultra- local in quantum mechanics (i.e., one-dimensional field theory with a coordinate A) 
and can be simply dealt with by an additive renormalization of mass m, as in the Euclidean space considered 
in [H]. 
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Note that we need to multiply ( 16. 8 j) by 1/2 to obtain the propagator G(x, x') of a neutral 
particle (i.e., particle = anti-particle): 



1 1 POO 

G(x,x') = ^A(x,x') = ± J dTe"^ T ^ ( m2 - i£ )- £ '/( 2T ) K(x, x'; T) 



(6.14) 



G(x,x') = i U/-g(y)d d yG(x,y) G(y,x'), (6.15) 



Since the propagator G(x,x') can formally be written as G(x,x') = i (x\ (□ — m 2 + 
ie)^ 1 \x') , one can easily show that G(x,x') satisfies the following composition law [9]: 

_d_ 

dm 2 

which is consistent with the asymptotic form of G(x, x') for large timelike separation with 
large mass: 

G{x,x') ~ e- imL(x ' x ' ] , (6.16) 

where L(x,x') is the timelike geodesic distance between x and x' . 

As a simple example, we consider a neutral particle propagating in a <i-dimensional 
Minkowski space with the metric 

ds 2 = -dt 2 + dx 2 . (6.17) 
Then the heat kernel can be calculated with the momentum representation as 



K(x,x'; T) = J 



dp 



(27l) d 

1 \d/2 r (x-x') 2 l 

Here, the first i reflects the fact that the Gaussian integral over p has the opposite sign 
of quadratic term to that for the other variables pi (i = 1, . . . , d — 1). Substituting this to 
( 16.141) . we obtain the following expression: 

-(d-2)/2 poo 

(47r) d / 2 

with z 2 = e' — i(x — x') 2 and a = i(m 2 — ie) . This integration can be easily performed, and 
we obtain 

j-(«J-2)/2 d _ 2 d _ 2 
X ') = ( 47r) d/2 2 * (V^ Z ) _ ~ K &=* (v^ *) 
rnf 1-2 d— 2 



g(g,a: / ) = * . ... / dss~ie-^- as (6.19) 



J/2 - (mVcr + ie') 2 (mVa + ie') , (6.20) 

where a = (x — x') 2 . This certainly agrees with the propagator (13. 29|) of a real scalar field. 
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6.2. de Sitter case 





L88| 


18 





In this subsection, we check that the in-out propagators (I4.5ip . (I4.87p . and 



Poincare patch 

We start from the following integral representation of the associated Legendre functions (a 
proof is given in Appendix [F]) : 



e **? /°° dA r (¥ + ^)r(^-iA) v; x % 2 w 

<W«)-* ^ dA r(iA)r( _ iA) L rfeZ > Re^>0j. 

(6.21) 

Since v 2 e has a positive imaginary part, we have 

j_ 9 d-2 d-2 

i f 00 f°° XT sinhfyrA) <*=i , .r,a. i3 , 

' dT / dA T ,<i . ' C.'^i^e 1 * ^ +A >. (6.22) 



2 2 ~f 7o io v/vr cos[tt(| - iA)] 
Here, u = Z — iO, and to obtain the second line we have used the identity 

-^2 d-a , 9 N ^a l)T(iA- ^) *=! 

p. 2 ( u ) = T— (u - 1) 4 \ - / / C. 2 d l (u) . 6.23 

iA -i v ; v ; r(rf/2)r(iA + ±±) iA -V v ' v ; 

We thus find that the in-out propagator (14.511) in the Poincare patch has the heat kernel 
representation of the form 



1 



.3 



G(x, x') = - I dT e- 1 ^ 1 K(x, x'; T) , (6.24) 



K(x,x'; T) 



o 



{2Tl) d / 2 

r(^i + iA) r(^i - iA) ^ ;t (,2, (d -i Y 2\ 

r(iA)r(-iA) 7 iA-i/ 2 We 



o 

6 2 



£(S) f dA ^) af^Me^^). (6,5) 



18 See [15] (also [9]) for the direct evaluation of the random walk in dc Sitter space, which is based on the 
heat kernel for Euclidean AdS space obtained in [16] . 
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The above heat kernel certainly satisfies Eqs. ( I6.12p and (16.131) . In fact, Eq. (I6.12p can 
be shown in the following way: 





i — K(x,x'; T) 
e -i^ r(¥) 

. d+l 

47T 2 



cos[7r(f-iA)] 



■- [-(Z 2 - 1) d 2 z K(x,x'; T) -dZd z K(x,x'; T 
■~O x K(x,x'; T) , 



where we have used the Gegenbauer differential equation 

'd- l\2i 



(6.26) 



(1 - u 2 ) d 2 u C.^{u) - dud u C.^(u) 



a A 2 _«W = 0, (6.27) 



and the fact that the Klein-Gordon operator for functions of the de Sitter invariant f{Z) 
can be written as 



Df(Z) = (1 - Z 2 ) d 2 z f(Z) -dZd z f(Z) . 



(6.28) 



The initial condition (I6.13P can be shown to hold by using the heat kernel equation and the 
equality 



D x -m 2 + ie)G' m/out {x, x') = — ^= 5 d {x - x') 



(6.29) 



as follows: 



:5 d (x-x') = - / dT e - i ^ T (n x -m 2 + ie)K(x,x';T) 



-<J 



o 



iK(x, x; 0) . 



(6.30) 



global patch 

In a similar way, using Eq. ( 16 .22 j) , we can show that the in-out propagators (I4.87P and (I4.88P 
has the heat kernel representation of the form 

1 r°° 2 • 

G i ; l Z\{x,x')= l - dTe-^ T K i ^ t (x,x';T) (6.31) 

1~ even J Z Jq L even * 
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with 



"dA A [CZ^iM ~ C£±±(u-)] e 1 * ( A2+ (^) 2 ) , (6.32) 



(-i)^r(^) 

47T 2 sm(7ri/J </0 

^in/out/ /. 
- fi even 1 - 1 / 

= , i 2 / / dA A tanh(vrA) [Cj^ (u + ) + Cj^ («_) e 1 * ( a2+ ^ 2 ) . 

4vr— cos(ttz/) 7o ~ ~ 

(6.33) 



Equations (I6.12p and (16. 13[) also hold for these heat kernels, as can be shown in the same 
way as above. 



6.3. Relation to the Green function in Euclidean AdS space 



As has been pointed out in [To], the in-out propagator in the Poincare patch is directly 
related to the Green function in Euclidean AdS space through an analytic continuation. In 
this subsection, we demonstrate this equivalence with a precise numerical constant. 

d- dimensional Euclidean AdS space (EAdS^) is defined as the hypersurface in a (d + 1)- 
dimensional Minkowski space with the relation 

VMN Y m Y n = - £' 2 (M, N = 0, . . . , d) , (6.34) 

where £' is called the AdS radius. EAdS^ has two connected components. A frequently 
used coordinate system which covers only a single connected component is the Poincare 
coordinates (z, y l ) (i — 1, . . . , d — 1) that are defined by the following embedding: 

y° = *° + * 2 2 + \v\\ r = r^-^-'"' 2 . ( 6 .35, 

We here have chosen the component with z > 0. The metric then takes the form 

ds2 = ^ dz* + dy ■ dy 

z l 
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The Green function in EAdS space is known to have the following form (see, e.g., [TT])r 9 

-i7r(d-2)/2 

where Z'(y,y') is the invariant of EAdS space 



GWlM/O = ( ^, d/2 (Z' 2 (y,y') - Q^/ 2 (Z'(y,y')) , (6.37) 



Z'(y, y') = - £'~ 2 VMN Y M (y) Y N (y') = 1 + { " Z ' f ^ ^ , (6.38) 



and 



u' = ^ (^y^) 2 + m 2 f 2 . (6.39) 

Note that Z' is always larger than unity. 

The coordinate system (z, y l ) is related to the Poincare coordinates (77, x l ) of <i-dimensional 
de Sitter space through the analytic continuation 

z = e iE ^(-ri), y' = x\ £' = e iE ^£, (6.40) 

or equivalently, 

Y° = iX d , Y { =X\ Y d =iX°. (6.41) 

In fact, one can easily show that the metrics of EAdS^ and dS^ transform to each other. 
One also finds the relations 



V < = J(- +m 2£>2 = Mt IV _ m 2 £ 2 +i Q 



2 

= u e , (6.42) 
= z 2 + z' 2 + \y-yf = (- V ) 2 + (- V ') 2 -\x-x'\ 2 _ . Q 

2 Z Z' 27]7]' 

= Z(x,x') - iO, (6.43) 
with which the Green function on EAdS space can be rewritten as 



e 



-i7r(d-2)/2 



G E A dS (y, V) = - , n W2 (u 2 - I)"" Qi d ~% 2 (u) (u = Z(x, x') - iO) . (6.44) 



(2vr) rf / 



n 



19 In fact, solving the Klein-Gordon equation with a delta function source using the EAdS invariant 
Z' = Z'(y,y') , we see that the Green function is a linear combination of (Z 12 — l) - ^ -2 )/ 4 ^?"^ 2 ^') and 
(Z' 2 — l) _ ' d_2 ^ 4 Q^I^2 2 (^') ■ By requiring that the Green function damps at large separation (cluster 
property), only the latter solution is selected, as can be seen from the asymptotic forms of the associated 
Legendre functions [see (|E. 19|) and (|E.20|) ] . The normalization is then determined by requiring that the Green 
function coincides with that in Euclidean space for infinitesimal separation of y and y' [or Z'(y,y') — > 1] . 
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This coincides with the in-out propagator (I4.5ip in the Poincare patch of de Sitter space up 
to a phase: 

G^(x, x') = e -j^- {u 2 - 1)-^ Q^/f W (« = Z(x, x') - iO) . (6.45) 

As pointed out in [5], the Green function of EAdS space has no direct relation with 
the in-in propagators associated with a vacuum of de Sitter space for any a. We see that 
it is the in-out propagator (in the Poincare patch) which is actually related to the Green 
function of EAdS space. We thus expect that we can obtain a deep insight on the dS/CFT 
correspondence [TS] by analytically continuing the Euclidean AdS/CFT correspondence and 
by interpreting the result in terms of the in-out propagators (not of the in-in propagators). 

7. Discussions and conclusion 

In this paper, we have considered quantum theory of a free scalar field in a nonstatic space- 
time. We first developed a framework to treat a harmonic oscillator with time-dependent 
parameters, and then applied it to investigate a free scalar field in de Sitter space, both in 
the Poincare and the global patches. 

We have introduced the vacuum state at each moment ti as the instantaneous ground 
state of the Hamiltonian at the moment. We developed a calculation method to obtain the 
wave function tp(t; ti) associated with the vacuum. The in-in and in-out propagators are 
then obtained from the wave functions by sending the initial and final times to the past and 
future infinities. 

A major advantage of our prescription in defining the vacuum is that we do not need to 
introduce "positive-energy wave functions" that cannot be defined in a definite way for a 
spacetime with no asymptotic timelike Killing vector. 

We have applied our method to calculate the in-in and in-out propagators in de Sitter 
space. The obtained propagators take de Sitter invariant forms, and are consistent with the 
results known in the literature. What actually happens is that, when the time tj is sent to 
the temporal boundary, our wave functions ip(t; tj) may have different forms from those of 
the in- or out-wave functions in the literature, but the obtained propagators yet coincide 
with each other (see footnote dH of subsec. 14 . 3 p . 

As a new result, we have found that a finite massless limit exists for the in-out propagator 
in the Poincare patch (with minimal coupling). This is in contrast to the in-in propagator, 
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where the no-go theorem is known that no massless limit exists for the in-in propagators 
without breaking the de Sitter invariance [3J. The same functional form had been obtained 
for the in-out propagators without precise numerical coefficients in 0, US] from other ap- 
proaches, and the massless case also had been considered in [9] . However, one cannot discuss 
the existence of a finite massless limit without knowing the precise numerical coefficients. 
Indeed, our in-out propagator in the global patch diverges in the massless limit just because 
the numerical coefficient diverges. 

We have argued that our in-out propagators should coincide with the Feynman propa- 
gators obtained by the path integral with the ie prescription, and have given a numerical 
evidence for the coincidence. We have also shown that the in-out propagators in both the 
Poincare and global patches satisfy Polyakov's composition principle. These observations 
exhibit the naturalness of the in-out propagators. It should be interesting to study under 
which conditions the coincidence between the in-out propagators based on the instanta- 
neous vacuum and the propagators derived from the path-integral formalism is realized, 
and furthermore, to investigate whether the composition principle holds universally for such 
propagators. 

Our in-in propagator in the global patch has a finite value for m > (d — l)/2, but it 
diverges for m < (d — l)/2. It will be important to compare the in-in propagators with 
those obtained (numerically) by the path integral of the Schwinger-Keldysh type [HI H2] 
(see also [T9]). 

As an important application of our construction, we can consider the Unruh-DeWitt 
detector working in a finite time interval. In the literature, in order to study the thermal 
properties of de Sitter space, the Euclidean vacuum has been considered. Since it is a 
vacuum state at the infinite past in the Poincare patch, the detectors have to start at the 
infinite past, which is the situation hard to realize in actual experiments. On the other hand, 
since our method allows us to construct the a very far but finite past, one 

may investigate the inherent nonequilibrium property of de Sitter space, by taking account 
of the corrections to the in-in propagators coming from the effect of such finite-time start 
(which violates the de Sitter invariance) [20]. It would be also interesting to investigate 
some physical quantities such as the rate of vacuum decay at finite times on the basis of our 
formalism. 

As another future direction, it should be interesting to apply our method to quantum 
field theories in spacetimes with horizons, such as spacetimes with black holes and de Sitter 
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space in the static patch. For such a spacetime, one needs to carefully study the consistency 
of our formalism with boundary conditions at the horizon. It would be also interesting to 
consider the 0- vacua in the rotating Rindler-AdS spacetime [21]. Since the inequivalent 
spacetime foliations (parameterized by 0) cover the same portion of the AdS spacetime, 
one expects that the propagator is uniquely determined by performing the path integral. 
However, the in- and out- vacua still has a nontrivial dependence on the parameter labeling 
the Hamiltonians. Thus it would be possible to observe which is preferred by the path- 
integral formalism. 

It should be important to consider interacting fields in generic nonstatic spacetimes and 
to establish perturbation theory on the basis of our formalism. It will be also interesting 
to investigate the in-out propagators for gravitons, since our method can be applied to field 
theory of higher spins without any essential modifications. 



Acknowledgments 

This work was supported by the Grant-in-Aid for the Global COE program "The Next Gen- 
eration of Physics, Spun from Universality and Emergence" from the Ministry of Education, 
Culture, Sports, Science and Technology (MEXT) of Japan. This work was also supported 
by MEXT (Grant No. 23540304). 



A. Proof of Eq. (ESS) 



Setting t = T s in (I2.47P and using the hermiticity ql(T s ) = q s (T s ) and pl(T s ) = p s (T s ), we 
obtain 



(W p [f,g]) \W p [f*J] W p [f*,g]J \c 2 )~ \c 2 
Then, from f !2"^g|) and ( 12"32"j) . we obtain 

B ;w) = [crW^c(oft)-A( t ;oft , |. (A.2) 

a*(t)J \a(t')J \a(t' 

A straightforward calculation shows that 

A( f -0 = i| - W P^*( T -*)^( r -* / )] -W p [F(T s ,t),<p(T s ,t>)f 
W p [<p*(T s , t), cp(T s , t')} W p [ip*(T s , t), ^(T s , t')} 
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Note that detA(t;£') = — 1 due to the commutation relations a> (t)] = — 1 and 

[ct(t'), a(t')] = 1 (this can also be checked by a direct calculation). Then, setting t = t' = t 
in flA~2|) and f!A~3l) . we find that 



_ t W p [<p*{T a ;t ),<p{T a ;t )] t 



Wy^*(T s ;* ),</?Cr s ;*o)] U Wp[^*(T s ; t„), </?(T s ; t )\ 



a . 



(A.4) 



B. Asymptotically Minkowski space 

In this appendix, we investigate a well-studied case where spacetime is asymptotically 
Minkowski in both the remote past and the remote future. 



B.l. Setup 

For brevity we consider the two-dimensional spacetime with the metric 

ds 2 = a 2 (t) (-dt 2 + dx 2 ) , 

where the scale factor a 2 (t) now depends on time and takes the form 

2 „ 1 - tanht 2 1 + tanht 
a [t) = a h a\ . 



(B.l) 



(B.2) 



This spacetime is asymptotically Minkowski with scale ao in the remote past and with scale 
a\ in the remote future. By expanding a scalar field 4>(t,x) as 



(f)(t,x) = ^2^2 ^k,a(t)yk,a(x) 

as in section ffl the a^ionbeconj 



k>0 a 



S[<j>(t,x)]= jdtdXy/^jj --gTd^dv 



m — is 



k>0 a 

Thus, the correspondence with the ingredients of section [2] is given by 



q(t) = </> k , a (t) , p(t) = l, u(t) = \lujl 



1 — tanh t 



+ u> 2 



1 + tanh t 



(B.3) 



EE f^\[Ka{t)-(k 2 + rn 2 a\t))<g^)]- ( B - 4 ) 



(B.5) 
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Since ie plays no essential role in this appendix, we have eliminated it from the action. 
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where 

u = \jk 2 + m 2 ag , u>i = \jk 2 + m 2 a 2 . (B.6) 

We also introduce 

w± = - (wi ± w ) . (B.7) 

B.2. Wave functions 

The equation of motion takes the form 

= g + ^ + c, 2 g = g+(^i^ + ^i T ^tanht)g, (B.8) 

which can be solved analytically with the hypergeometric function. We set a pair of inde- 
pendent solutions {f(t),g(t}} as 

,\ _ £n iwi/2 /l -f £\ -iwo/2 / 



= (ht) ^(-^->i-^-;i + ^o;^), (b.io) 

where £ = tanht and F(a,b;c; z) is the hypergeometric function. Their asymptotic forms 
for t — to ~ — oo (or £ = £ ~ — 1) are easily found to be 

Zo = /(*o)~e-^\ <7 = <7(*o)~e iw °\ (B.ll) 

and the weighted Wronskian W p [f, g] = pit) W[f,g](t) is found to be 

W p \f t g]=2iu . (B.12) 

From this we find that the functions in (I2.49P have the asymptotic forms 

u ~0, u 2iu e- iuloto , (B.13) 

v ~ 2iw e iwoi °, w ~0. (B.14) 

The wave functions in the limit t — > — oo then take the form 

^ ^^to-i^t-i*- Iog(2oo B ht) F ^ iw _ j 1 + ; 1 _ iwo ; 1+t 2 allht j , (B.15) 

(p{t-t Q )~-^ =e - i »^g{t) 

^ e -i W oto+i W+ t+i W -iog(2co«ht) F / r _ i£i; _ > i _ iw _ ; i + iWo ; 1 + tanht ) . (B.16) 
2ujq V 2 / 
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In order to calculate the asymptotic forms of various functions for t ~ +00 , it is conve- 
nient to rewrite f(t) and g(t) by using the formula 

/ , n T(c)T(c — a — b) , . 
F(a, b] c; z) = - ) ' \ -( F(a, 6; a + b - c + 1; 1 - z) 



T(c- a)T(c-b) 
r(c) r(a + b - c) 



+ 



(1 - z) c - a ~ b F(c-a, c-b; c-a-b+1- 1- z). 



T(a)T(b) 



(B.17) 



We then obtain 



/(*) 



5/1 -C 



iojo/2 



F( 1 — icj + , — iu + ; 1 — iu\ 



ujo\ x / 2 / l -C \" icJl/2 / 1 + C \ 1 



- 



ioj /2 



a*F[ iu-, 1 + icj_ ; 1 + iui ; 



i-C 



aFl — ia>_, 1 — icj_ ; 1 — iui 



i-C 



i-C 



5*/!-C 



where 



a = 



uji\ 1/2 r(l + ujq)T(uji) 
o; / T(l + iw+)r(io; + ) 



F( 1 + iu;+, iw+ ; 1 + iwi 







i-C 



wa 1 / 2 T(l - iw )r(ia;i) 
r(l + iu-)T(iu-) 

With these, the asymptotic forms for t\ ~ +00 can be obtained easily as 

fx = f(ti) ~ f^) 1/2 (a*e- iwitl - ^e^ 1 ) , 



»i = y(ti)~(^) 1/2 (fie 1 ^-^ e -^) 



and 



Ml 



(^)' /2 (-2.M) e 



itiJiti 



Ml 



1/2 



-2ia*wi)e 



-iuitl 



iojifi 



^V /2 (2i / S*o; 1 )e- ia ' 1 * 1 



We then find that the Bogoliubov coefficients take the asymptotic forms 
a{h,to) ~ a e-^oto-uiti) = ai (a{t u t ) ~ af) , 



(B.18) 

) 

(B.19) 

(B.20) 

(B.21) 
(B.22) 

(B.23) 
(B.24) 

(B.25) 
(B.26) 



They coincide with the well-known values in the literature (see, e.g., in [7]) up to a phase. 
It is easy to see 

,^ sinh 2 (7ru; + ) sinh 2 (7ru;_ 



sinh(7ra;o) sinh(7ru;i) 



sinh(7ro;o) sinh(7ru;i) 



(B.27) 
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and thus the relation |ai| — =1 actually holds. 

Using the asymptotic forms of U\ and V\ , we can calculate the wave function (/?(£, t{) for 
t\ — > +00: 

1 



af(t) + 0g(t) 



This can be further rewritten by using Kummer's relation 

F(a, b; c; z) = (1 - z) c ~ a ~ b F(c - a, c - b; c; z) 



into the form 



p ia;iti-ia; + t~ia>- Iog(2cosht) F f • I + [ u _ ■ 1 



1CU 



1 — tanh £ 



(B.28) 



(B.29) 



(B.30) 



This certainly coincides up to a phase with the positive-energy wave function in the re- 
mote future given in [7J. One can easily see that ip(t;ti) actually has the form ip(t;ti) ~ 
(1/ y/2u)i) e -11 ^* - * 1 ) when t is also very large. 



C. Propagator in Minkowski space 

In order to evaluate the integral (I3.28p . we introduce the polar coordinates for the wave 
vector as 

dk 2 = dk 2 + k 2 (d6 2 + sin 2 9 dn 2 d _ 3 ) , (C.l) 

where 9 is chosen such that 9 = corresponds to the direction x — x' , i.e., k ■ (x — x') = 
k | x — x' | cos 9 . The volume element is then given by 

d^ 1 ^ = dk k d ~ 2 d9 sin d - 3 9 dfl d _ 3 = dk k d ~ 2 d cos 9 (1 - cos 2 9) {d - i)/2 dQ d _ 3 . (C.2) 

Equation f!3.28j) then becomes 

IO I f°° f l 
G(x,x') = d ~Z / dkk d - l G k (t,t') / ds(l-s 2 Y d - A ^ 2 cos(k\x-x'\s) , (C.3) 
{^) Jo J-l 

where is the area of the unit sphere in n-dimensional Euclidean space, = 

J df2„_i = 2 7r n / 2 /r(n/2) . The integration with respect to s = cos 9 can be carried out by 
using the formula (8.411-8 in |22j ) 

J ds{\- s 2 y cos{zs) = ^T{v+l)(^ U ~~ 2 J u+ i_{z) [Rev>-1], (C.4) 
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and we obtain 

x') = — / dfc fc~ Gk(t, t') Jd-3 (k\x — x'\) 

2 -(d+l)/2 n -(d-l)/2 roo k (d-l)/2 



\(d-3)/2 



/ dk e -i^ (*>-*<) J d _ s (k\x-x'\) . (C.5) 

Jo 2 v y 



Although the integration does not converge in the UV region, this expression still has a 
meaning as a distribution and is analytic in the lower half of the complex t> — t< plane. 
In order to make the integral manifestly converge, we introduce a small imaginary part by 
replacing t > —t < by t > — t< — iO which does not alter the analytic property. Thus, what we 
need to evaluate is the integral 

9-(d+l)/2 -(d-l)/2 roo Ud-l)/2 . 

G(x,x') = - -jT^r- / dfc- e- e7r/ 2 <^>-*<-°) J d . z {k\x - x'\) 

\x - x'\ (d - 3)/2 Jo —V 1 1; 



2 -(d+l)/2 ^-(d-lJ/2 
_ U -i(d— 3)/2 



/oo 
dA (A 2 - 1)^ e - mAe " /2 (*>-*<~ i0 ) (k\x -x'\), (C.6) 



where, assuming m > , we have set A = Uk/m = \/k 2 + m 2 /m to obtain the last expression. 
Then applying the formula (6.645-2 in 



1 

dA (A 2 - 1)2 e" QA J v (/3 V^^i) 

(a 2 + /3 2 )~ f ~^ AT I/+ i( v / a 2 + ^ 2 ) [Rea>0,/3eM] (C.7) 



with i/ = (d — 3)/2 , a = m e 17r ' 2 (t> — t< — iO) , /3 = m |cc — x'\ , we obtain 

m (d~2)/2 d _ 2 

G(x, x') = - . ... [e^ (At 2 - iO) + As 2 ] 4 Km (mJe^(At 2 - iO) + Aa; 2 ) , (C.8) 

(27Tj ' 2 

where At = t — t' and Aa; = x — x' . This expression can be further rewritten by separately 
investigating the cases for different sign of a = (x — x') 2 = —At 2 + Ax 2 . 

(1) spacelike (a > 0) : 



The modified Bessel function is readily evaluated as if( ( i_2)/2(?™\/e 17r (Ai 2 — ^0 + Aa; 2 ) 
K(d-2)/2( m V&) , and we have 



m (d-2)/2 _ d _ 2 



G(x,x') = ^y /2 o ^ ^(m^). (C.9) 
(2) timelike (a < 0) : 

By using the relations e i7r (At 2 - iO) + Aa; 2 = e' 17r (-a) and K u (e in/2 z) = - (in/2) H { y\z) , 
we have 

G(x, x>) = \ (2 ^ 7 L (-o-)-" # S (m V^) • (CIO) 
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(3) null (a -» 0+) : 
By using the expansion 



^(THi+oc* 2 )) 



2 V2 



we obtain 



(C.ll) 



(C.12) 



4^/2 

The right-hand side actually coincides with the massless propagator. 

It is easy to see that all of the expression for three cases can be derived from a single 
expression, 

m (d-2)/2 



G(x,x') 



(2n) d / 2 (a + iO) 



(d-2)/4 — 



Kd-2 (m\/ a + iO) . 



(C.13) 



D. Proofs of Eqs. ( 14471 ) and (ESI) 



In order to show Eq. (I4.47p . we start from the following equation (6.578-10 in |22j): 



dxx^ +l K u {ax) K v (bx) 



ex) 



^c»T{n + u + l)T{ii-u + l) ( 2 (m+ |) 



23/2 ^ 6)^+1 



(«»-!)■ 



^-1/2" 1 



a 2 + o 2 + c 2 f 

it = ^ , Re (a + 0) > |Imc| , Re (/1 ± z/) > — 1 , Re/x > — 1 

2a 6 



(D.l) 



By setting a = e 1?r / 2 a and b = e 17r / 2 6 with —n/2 < arga < ix and —n/2 < argfr < ix , and 
by using the identities 

2i 



01 



e-to/Jjf/*^ #(2)^ 

7T 

the equation is rewritten to the form 

dx x^ +1 (a x) {b x) (c x) 



x) = — e 

TV 



^ /2 K v (bx), (D.2) 



u 



_ a 2 _ 5 2 _)_ c 2 

— , Im (a - 6) > |Imc| , Re (jj, ± v) > — 1, Re/i>-l . (D.3) 
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We substitute to this x = k, a = —r] < — ie f , b = — r] > — ie, c = \x — x'\ and /x = (d— 3)/2 where 
the convergence of the integral requires the inequality e > e' . We then obtain Eq. (I4.47P 
with u = —Z(x, x') + i for infinitesimal e and e' . 

In order to show Eq. (I4.5ip . we use the following equation (6.578-11 in [22]): 



dx K v (a x) I v (b x) 



ex) 



o 



fi -i7r(/i+|) 1 x , 

r (n 2 -l)- ( ^ S5 /2 ( H ) 



cr e 



27r (dtV)' 1 - 



Q 2 _)_ 5 2 _|_ Q 2 A "I 

u= ; , Re a > |Refe| + |Imc| , Re/i>-l, Re(/i + z/)>-l . (D.4) 

2d b J 

By setting d = e 1?r//2 a and 6 = e 17r / 2 6 for —n/2 < arga < 7r and — n < argb < n/2 , and by 
using the identities 



2i 



H [2 \ax) = —e^ u/2 KJax) , J u (bx) = e~' 17ru/2 iJbx) 

TV 



(D.5) 



the equation is rewritten to the form 



dxx^ +1 H {2 \ax) J u (bx)J^(cx 
a 2 + b 2 - c 2 



e -2iTT(fl+±) 



( M 2 _ iy?in+k) , 



u 



2a b 



7 r 3 / 2 (aby+ 1 v " v ^ u - 1 ' 2 

Ima) > |Imfe| + |Imc| , Re/i>-l, Re(/i + ^)>-l . (D.6) 



We substitute to this x = k, a = —r] < — ie, b = — ?7> — ie', c = \x — x'\ and fi = {d — 
3)/2 , where the convergence of the integral requires the inequality e > e' . We then obtain 
Eq. (I4.5ip with u = Z(x, x') — i for infinitesimal e and e' . 



E. Associated Legendre functions and the addition for- 
mulae 

In this appendix, we give several formulae of the associated Legendre functions which are 
used in subsec. 14.31 For details of the associated Legendre functions, see [22] and [23] . 

The associated Legendre functions V%{z) and Q„(z) are defined over the complex z-plane 
other than the cut along the real axis to the left of the point z = 1 (running from — oo to 1), 
while the associated Legendre functions P^(x) and Q„{x) are defined only on the interval 
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-1 < x < 1: 



KM 

Pt{x) 



z + 1 1 1-2 

i< -z/, v + 1; 1 - /i; — - — 



r(i -n) \z - i 

2"+ 1 r(i/ + 3/2) 



(E.l) 



(* 2 -i)*f( 



z/ + /i + 2 Z/ + /U+1 



3 1 

2' 1^ 



e - 2 iwp*( x + io) + e"3 i7r ^(x - iO) 

i n + x\% 



F[ —u, v + 1; 1 - //; 



1 -x 



r(i - /i) vi -x 

1 Q£(x + iO) + et i7r " Q£(x - iO) 



-17T/i 



e 2 



7T 



2 sin 7r/i 



cos7r M P^)- ^ + ^j| P^(x) 



(E.2) 



(E.3) 



(E.4) 



functional relations 

The four functions V(z), Q(z), P(x) and Q(x) are related to each other as (3.4 and 3.3.1 in 



i7r c 

e 2 



-lTTfl 



V?(x + iO) = V^x - iO) = P£(x 
Q!t{x + iO) ± e 1 ^ Q»(x - iO) 



2Q^ 



-ivrP^OO 



es(-*) 



-e^Q^) [ImoO] 
- e-' 11fV Q£(z) [Imz < 0] 



(E.5) 
(E.6) 

(E.7) 



We also have (8.73 in [22] ) 

r^-zx + i) 



P^(x) 



r(z/ + /i + 1) 



7>£(z) - -e"^ sinvr/i Q£(z) 

7T 



r(i/ + ^ + i) 



r(z/-/i + 1) 



-x 



Q' 



X 



X 



r(i/ + // + i) 

I>-/i + l) r 

r(z/ + /i + 1) 
r(i/-// + 1) 
r(i/ + /x + i) 

l 

sin 7r(i> — ji) 



COS7T/iP^(x) sin7r/iQ^(x) 



7T 



Sin 7T/X P^(x) + COS ir/i (x) 



cos7rz/P^(x) shi7rz/Q^(x) 

7T 



[ — 7r cos7r^ Py(x) + sin n(is + ii) Q^(x 



(E.8) 
(E.9) 
(E.10) 
(E.ll) 

(E.12) 

(E.13) 
(E.14) 
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Their Wronskians have the forms (8.741 in 



d.r 



i r> + /x + 1) 



P^(x) P£(x) 



i - x 2 r(f - fi + 1) 

1 2 sin 7r/i 



1 -x 2 



7T 



We also have (8.733-1 in [22] 



(1 " T" P 2C=) = + P£(*) - (" - A* + 1) P£h(*) 



dx 



(1-x 2 ) — Q£(x) = (z/ + l)xQ£(x) 



/i + l)Qe +1 (x) 



(E.15) 
(E.16) 



(E.17) 
(E.18) 



asymptotic forms 

The associated Legendre functions have the following asymptotic forms near boundaries 
(3.9.2 in [23]): 



i r ^ + i/2 

, i r(-v- 1/2) 
2-u-i n - 2 y /) z -u- 1 

r(-v-fi) 

i 

7T2 — - : — - Z 



[Rev > -1/2] 
[Reu< -1/2] 



i7T/J. 0-J/-1 „i T(Z/ + /i + 1) 

~ e ^ 2 



Q£(x) 2"" 1 cos 



r(z/ + 3/2) 
(l-x 2 ) _i , 

{jxv)Y{y) (l -x 2 )~ f [Rei/>0] 

r(fc + i/ + i) 



(E.19) 

(E.20) 
(E.21) 
(E.22) 



V 1 T{v+l)T{k-v + l) 
2 V sin(vrA;) T(v 



(l-x 2 ) 2 [fcez; 



7T 



(1-x 2 )" 



2"" 1 cos(ttA;) r(z/) (1 -x 2 ) 2 

2-^ 1 7rr(A; + z/ + l) 
sin(vrA;)r(z/ + l)r(A;-z/ + l) 



[jkeZ+1/2, Rei/ > 0] 

(E.23) 

[ Jfe 6 Z , Re i/ > ] 



(l-x 2 ) 2 [jfeeZ + 1/2] 



. (E.24) 



addition formulae 

We find the following formulae, which are useful in obtaining the propagators in the global 
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patch: 



17T(COS</?i COS<y2 2 J 

2(d- 2) sm(iriy] 



-i oo 



J2( 2L + d - 2 ) P fc "( sin Vi) P fc( sin ^2) C L 2 (cos 1 



L=0 



;„ d-2 

- e" 17r — 



2 (27r) 2 sin(7r^) 



« - 1)-^ ezi(«+) - («- - 1)-^- ezi(«- 



(d : odd) 



(E.25) 



. / \ d—l OO 

l(COS(^i cos<£> 2 j 2 



(d — 2) COs(7TZ/) 



^(2L + d - 2) "(sin </?i) Q£(sin <^ 2 ) C L 2 (cos 0) 



L=0 



ie~ ln — 



2 (27r) 2 C0S(7TZ/) 



(d : even) 



(E.26) 



. / \ d-l OO 

i(cosv?i cosy? 2j ) 2 



(d-2)\n d . 



;_ d-2 
le -17T — 



J^(2L + d-2) 



rCfc — 1/ + 1) d-2 

J P^( S m^)Q^(sm^)C^ (cosfl) 



2(2tt 



L=0 

-i7Tl/ / ,2 iN-^M^^ \ 1 Attu /.2 



QZ-i («+) + ei ™ - 1)"^ si_i(«- 

(d : even) , (E.27) 



_ . / , d-l OO 

21(008^ cos^ 2 ) 2_ ^ (2L + d _ 2)Qr(sin ^ )Q , (sin ^^ (cos( 



ir(d-2)\n d _ 



L=0 



:_ rf-2 



2(2*; 



+ 



1 7T 



COS 7TZA 



[e" 1 ™ {u\ - 1)-^ + {ul - 1)-^ 7^(u_ 



2(2tt) 



7T 



COS 7Tl> 



« - 1)-^- Q„\(U+) - {ul - 1)"T- Q\{U- 



[« - p„ \(-u + ) - (u 2 - i)-"?!iK 



(d : even) , (E.28) 



where —ir/2 < ip 2 < (fi < tt/2 , < 9 < ix , k = L + (d - 3)/2 , and 

u±(tpi,tp2,0) = -^(v?i,^2,6 | ) ± iO with Z(tpi,tp 2 ,9) = sm ^1 sm ^2 + cos 6> ^ (E.29) 

COS<£>i COS<£> 2 

We prove Eqs. (1E.25|) and (1E.26[) for the rest of this appendix. Equation ( 1E.27|) can be 
proved in a similar way, and ( IE. 281) is readily obtained from ( IE. 261) and ( IE. 271) . 
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We start from Eq. (12) of [23] 



(sinh 7) 2 Q \ (cosh 7) 



Hr(Y)e"h") (sinh ft sinh/3 2 )V 

00 

^(2L + c/-2)p- 1/ (cosh/3 1 ) Q£(cosh/3 2 ) C^(cos0) , (E.30) 



L=0 



where 



, „ „ „ N cosh /3i cosh /3 2 — cos a 
cosh7(/3 1 ,/3 2 ,0) = . , , Re/3 2 > |Re /3 X | + |Im0| . 



(E.31) 



sinh Pi sinh /3 2 

Both sides of Eq. ( 1E.30I) should be understood as the quantities that are continued analyt- 
ically from the region where Pi, P 2 , and 9 take all real values (for which cosh 7 > 1). We 
reparametrize the variables in Eq. flE.30j) as 



ft 

7l 



±i 



7T 



2 < if 2 < <fl < , 



(E.32) 



(E.33) 



(E.34) 



„J+e Q (a = 1,2) 
< £1 < e 2 < 1 

and only keep the contributions from e a to the linear order. We then have 
cosh P^ = sin ip a ± \e a cos <p a , sinh P^ = ±i cos cp a + e a sin <f) a 

and 

cosh7 ± = cosh7(/3 1 T ,/3 2 ± ,6') = -Z(ipi, (p 2 , 0) + iO(ei,e 2 ) . 

If we fix the parameters £1 and e 2 , and vary <£> 2 , and 9 within the regions —tt/2 < y? 2 < 
<f 1 < 7r/2 and < 9 < tt , then cosli7± ranges in the region depicted in Fig. |6j 

In the following, we divide the parameter region of tpx, <p 2 and 9 into three parts, where 
Z takes values in (1) Z > 1, (2) 1 > Z > — 1, and (3) Z < —1, respectively. We then derive 
a simpler expression of Eq. ( 1E.30|) for each case, and show that the obtained expressions for 
the three cases can be summarized in the form flE.25j) and (IE.26j) . 

(1) z>\ 

In this case, cosli7± = — Z ± iO as can be seen from Fig. I6l 22 l Thus, we have 

sinh 7± = e ±iw (Z 2 - l) * , (E.35) 



Q „l i(cosh7 ± ) 



3 =Fi7r(i/-i) 



(E.36) 



21 





22 



Equation (|E.27[) can be proved by replacing (cosh ^2) in (|E.30|) by Q k ^ (cosh ^2) with the help of 
I 

£\ and £2 are to be taken to zero with keeping e 2 > £j.. 
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l.o r 



0.5 - 



cosh 7^ 



Figure 6: Schematic view illustrating the range of cosli7 + for E\ = 0.001 and 
62 = 0.02. Here, tpi and ip 2 run over the range — 7r/2 < tp 2 < (pi < 7r/2, and 
9 runs over its full range < 9 < n . The range of cosh 7_ can be obtained by 
turning the figure by 180° over the horizontal axis. 



where we have used Eq. (IE. 71) and the fact that Q„{z) does not have a cut in the region 
Rez > 1. Then, Eq. ( IE. 301) becomes 



d-2 d-2 



oo 

x J^(2L + d - 2) "(sin (f!) Q£(sin (p 2 ± 10) (cos 9) 



(COS(/?i COSy? 2 ) 2 



(E.37) 



L=0 



By taking the difference between the above equations with the upper and the lower signs, 
we obtain 



2i sin 



r /d-3 




7T 


-)] 


L V 2 





d-2 d-2 



d / dt — 2\ ■ d-2 

i7r22~ 2 r 



{Z 2 -iy 4 Q u \(Z) 

; d-2 

e 2 (cos</?i cos(/?2) 



d-1 
2 



X 



J2( 2L + d-2) Pfc "(sin pi) Pfc(sin <p 2 ) C L 2 (cos 0) , 



(E.38) 



L=0 
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where Eq. (1E.6j) has been used. Similarly, by taking their sum, we obtain 



2 cos 



2~ 1 r 



" (d — 3 




K 2 - 









d-2 



(Z 2 -l)~ 4 ^\(Z) 



e 2 (cosy?i cos y? 2 ) 



d-l 
2 



5^(21 + d - 2) P" ^(sin ^) Q£(sin <p 2 ) c/ (cos i 



L=0 



The right hand side of Eqs. ( IE. 251) and ( IE. 261) can then be written as 

'd-2 



iSin7r( ^ } (Z 2 - 1)-^ Q^(Z) (d : odd) 

V 2 



(2tt)I 



1 COS7T( tan7ri/ 



d-2 d-2 



(27T 



e^ i7r — (Z 2 -l) 4 Q^i(Z) (d:even). 



(2) 1 > Z> -1 



(E.39) 



(E.40) 



In this be seen from Fig. [6l cosh 7+ crosses the branch cut between — 1 < 

z < 1 from above and move to another Riemann sheet. On the other hand, cosli7_ crosses 
the branch cut between — 1 < z < 1 from below. Thus, in this region, we have 



sinh 7± = e ±i 5 (l - Z 2 ) 3 , 
Q^(cosh7 ± ) = Q~f k (-Z± iO) = e^ (1± 5) 



(E.41) 



d-2 177- d-2 

Q u \(-Z)T-P u \(-Z) , (E.42) 

^2 z 2 



where we have used Eq. (IE.6p . Then, Eq. (1E.3QI) becomes 



;_ d-2 



22-^ r 



(1-z 2 )- 

d-2 



d-2 
4 



d-2 i 77- d-2 

Q^-z^-Pj^-z) 

^2 z 2 



. l "(-"T8+ 4 f 



) 



COS(/?i COS(/?2j 



J^(2L + d - 2) P^(sin^ 1 ) Q£(sinp 2 ± iO) c/ {cos9) . (E.43) 



L=0 



By taking the difference of the equations with the upper and the lower signs, we obtain 



d-2 d-2 



(i-z 2 ) ~p \(-z) = 22- 2 r 



d-2 



(COS tfii COS</9 2 ) 2 



X 



J2( 2L + d-2) Pfc"(sin^ 1 ) Pfc(sin^ 2 ) C^ 2 (cos#) . (E.44) 



L=0 



Similarly, by taking their sum, we obtain 

d-2 d-2 rf _ /rf_2 



(1 - z 2 )"^ Q v Ti(-^) = 2l " 2 r (^) ( cos ^ cos ^) 



d-l 
2 



X 



J^(2L + d - 2) Pfc ^(sin </?i) P£(sin tp 2 ) c/ (cos 0) . (E.45) 



L=0 
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The right hand side of Eqs. (1E.25[) and ( 1E.26j) can then be written as follows: 

\Tl , d — 2 d — 2 



— (l-Z 2 ) 4 P \(-Z) (tZ: odd), 
2( 27 r)fsinM ^ (E46) 

^ ; : (1 " 4 Q^i(-^) (d ■ even) . 
(27r) 2 cos^z/J 2 



(3) Z< -1 



In this case, from Fig. [61 we know that, in the region Re 2 > 1, cosli7 + runs above the 
real axis, or below the real axis in the next Riemann sheet after passing through the cut on 
— 1 < z < 1 from above. On the other hand, in the region Kez > 1, cosh7_ runs below the 
real axis, or above the real axis in another sheet after passing through the cut on — 1 < z < 1 
from below. 

When cosh 7-|- = — Z ± iO, we have 

sinh 7± = (Z 2 - 1) 1 , (cosh 7± ) = (-Z) , (E.47) 



and then the following equation is obtained: 

d-2 d-2 



{z 2 -iy 4 q v \{-z) 



22- 2 r 



d-2 



2 

^(2L + d -2) Pfc^sin^x) Qfc(sin^ 2 ± iO) C^cosfl) . (E.48) 



oo 
L=0 

By taking the difference and the sum of the above equations with the upper and the lower 
signs, we obtain 

d 2 /d — 2\ in d-2 , ,d_i 



= - i7r 2 2 Z T[ — — ) e 17r 2 (cos <p x cos <p 2 ) 2 



oo 

x J^(2L + d- 2) P- V (sin^ 1 ) P£(sin<^ 2 ) C^(cos0) , (E.49) 

L=0 

/ 9 \-— — / x d 1 /d — 2\ . d-2 . , d-i 

2(Z 2 -l) 4 Q^i(-Z) = 22~ 1 r(^^— Je 1 "— (cos^! cos^ 2 ) — 

°° d-2 

x ^(2L + d-2)P^ v (smtp 1 )Q v k (smtp 2 )C^~(cos9). (E.50) 

L=0 

On the other hand, when cosh 7 -|- = — Z =p iO, we need to evaluate the functions sinh 7-1- and 

d-2 

Q %(cosh7±) on a new Riemann sheet, since cosh7± has already crossed the branch cut. 

u ~ 2 

We then have 

sinh 7± = e ±i7r (Z 2 - l)* , (E.51) 

Q i X 1 (cosh 7± ) = e ±i7r ^ QZd-Z) =F iTre^T^-Z) , (E.52) 
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and Eq. (1E.30|) takes the form 



d-2 r d — 2 j_„ , d-2 

(Z 2 - QJ\(-Z) T i7re ±i ^ (lTl) V^\(-Z) 



d n „ / d — 2 

22-2 r' 



oo 

J2( 2L + rf - 2) Pfc"(siii^i) Qfc(siri(^ 2 ± iO) C^(cos0) . (E.53) 



L=0 



By taking the difference of the above equations with the upper and the lower signs, we 
obtain 



d-2 



e i7r V(cos<^i cos^ 2 ) V 5^( 2L + d ~ 2 ) P r( sin< ^i)Pfc(sin^)Ci 22 (cos0) . 



L=0 



(E.54) 



In particular, in odd dimensions, we have 

i 7t — — — / 

e 2 (cOSy?i COS<y2 2 )' 



o = 22 _2 r 



x 



J^(2L + d - 2) P*T "(sin <px) P£(sin y> 3 ) Q 2 ' (cos i 



(E.55) 



L=0 



which is equivalent to Eq. ( ]E.49j) . Similarly, by taking the sum of Eq. (1E.53|) with the upper 
and the lower signs, we have 

(z»-ipr 

'd — 2\ -^d- 



d-2 j_„ / d — 9 \ d-2 

\ (-Z) - Tre 1 ^ sin (—tt) P % (-Z) 



2t-T 



e iw T a (cos^ 1 cos(p 2 ) V 5^(2^ + d - 2)P^(sin^i)Q£(sin^ 2 )C , i 2 (cosfl) . 

(E.56) 



L=0 



In particular, in even dimensions, we have 

d-2 d-2 



22- 2 r 



d-2 



e 2 (^cos^i cosy2 2 j 2 



J^(2L + d-2) P k "(sin </?i) Q^(sin </? 2 ) C L 2 (cos#) 



L=0 



(E.57) 



which is equivalent to Eq. (1E.50|) . Thus, when Z < — 1, the right hand sides of Eqs. (1E.25[) 
and ( 1E.26I) always take the forms 



(d : odd) 



(2ir) 2 cos(7rz/) 



;_ d-2 



d-2 d-2 



[Z 2 -\Y * Q\(-Z) (d: 



even 



(E.58) 
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We thus have obtained simplified expressions for Eq. (1E.30h for three different regions of 
Z in the form Eqs. ( IE .40 p . (|E.46j) . and (IE.58[) . One can readily see that three equations can 
be obtained from Eqs. ( IE.25j) and (IE.26j) . This completes the proof of our assertion. 



F. Integral representation of the associated Legendre 
functions 



In this appendix, we give a proof of Eq. (I6.2ip [we write it again here for convenience 

r(4=i + ia) r(fci - iA) v; x % 2 



d-2 
2 



dX 



r(iA)r(-iA) 



z/ 2 + A 5 



[deZ, Re v > 0] , 
(P.l) 



treating the odd and even dimensional cases separately. Our discussion is heavily based on 
the derivation of the heat kernel in Euclidean AdS space performed in [16] (see also |15j). 



odd dimensions 

When d is odd, using the direct relation between the associated Legendre functions and the 
Gegenbauer functions, 

2%V/ 2 r(^i) 



v 



w- 



[U 



sin^-iA^r^ + iAK^-iA) 



1)—C 



d-i w , 



we can rewrite the integral on the right-hand side of (IF.ip as 



u 



°° HA r(^ + iA)r(^-iA) <y 

T(iA)r(-iA) ^ 2 + A 2 



A i=i 



By further using the identities for the Gegenbauer function 

2 1_n r d 

iACT A _ n (cosh 7 ) 



T(n) 



d(cosh7) 



cos(A7) 



with nonnegative integers n , Eq. (1F.3I) can be rewritten as follows: 



2 l/2 e i7r 



d-l 



(cosh 2 7 - l) 4 
(cosh 2 7-1)^ 



7rV2 



-d(cosh7). 
d 

d(cosh7) 



^ /' 00 dA cos(A7) 



d-l 
2 



z/ 2 + A 2 
T7 ^ [Re 7 ^0, Rez/>0 



Q^* (cosh 7) [Rei/>0]. 



(F.2) 



(F.3) 



(F.4) 



(F.5) 
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Here, in order to show the second equality, we have used the formula [3.723-2] of 

cos(A7) 



7T 



dA 



— e^ u [Re 7 >0, Reu>0] 



v 2 + A 2 



7T 

27 



e lu [Re7<0, Re^>0] 



(F.6) 



and to show the third equality we have used the following formula [see (A. 20) of 

d 



Qu-i/2 ( cosh 7) = < 



2V2z/ 
7r-'" e i,r 2 



(cosh 2 7 - l) 4 



1/2 giTT^ 

(cosh 2 7 -l) 4 



-d(cosh7) 
d 



d-l 
2 



-7 V 



d(cosh7) 



d-l 
2 



3 7 " 



[Re 7 > 0] 

[Re 7 < 0] 

(F.7) 



even dimensions 

When d is even, the equality 



rfiA- — + l) d-2 

iA -5 v ; r (iA + ^i) iA -2 ; 



holds, and we can show fIF.ip as follows: 

r(^i + iA)r(^-iA) v; x % 2 (u) 



dA 



r(iA)r(-iA) 



z/ 2 + A 2 



dA 



A tanh(7rA) 



V 



+ ^2 ' U-l/2 



«/d\ 5 r f,J tanh(TrA) 



dA 



(-l)T-(tl 2 -l)~^ ^ 

(-l)^^-!)^ (A) 

■ d-2 ±ll 

e"" r " r ^i/ 1 («). 



+ A 2 



P 



1A-1/2W 



(F.8) 



(F.9) 



where we have used the identities for the associated Legendre functions with integer order, 



Jn 



/2 



du 

and the formula [7.213] of [22]: 

A tanh(7rA) 
^ 2_ 



du r 



Q„i(u) = Q n v _ l (u), (F.10) 



dA 



PiA-i/ 2 («) = Qu-i/2(u) [Re i/ > 0] . 



(F.ll) 
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